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Abstract 



Consider Glauber dynamics for the Ising model on the hypercubic 
lattice with a positive magnetic field. Starting from the minus configu- 
ration, the system initially settles into a metastable state with negative 
magnetization. Slowly the system relaxes to a stable state with positive 
magnetization. Schonmann and Shiosman showed that in the two dimen- 
sional case the relaxation time is a simple function of the energy required 
to create a critical Wulff droplet. 

The dilute Ising model is obtained from the regular Ising model by 
deleting a fraction of the edges of the underlying graph. In this paper we 
show that even an arbitrarily small dilution can dramatically reduce the 
relaxation time. This is because of a catalyst effect — rare regions of high 
dilution speed up the transition from mirms phase to plus phase. 
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1 Introduction 

1.1 Metastability in the Ising model 

Consider Glauber dynamics for the supercritical Ising model on the hypercubic 
lattice {d ^ 2) started in the minus configuration but with a positive external 
magnetic field h. Aizenman and Lebowitz predicted that the model initially 
settles in a metastable minus phase, eventually relaxing to the plus phase on a 
time scale that grows exponentially with l/h'^~^ [T]. 

To be more precise, let /3 denote the inverse-temperature and let /3c denote 
the critical inverse-temperature. Suppose j3 > (3^- Let ix^^fjT denote the plus 
and minus phases of the equilibrium Ising model. Start the Glauber dynamics at 
time with all vertices initially taking minus spin. Let (t^'~ denote the state of 
the Glauber dynamics at time t. With (3 fixed, let ft- —> with t — exp(A/ft'*~^). 
A heuristic argument suggests that if Ai is sufficiently small and A2 is sufficiently 
large then for every local observable /: 

(i) E[/(a°'-)] if A<Ai. 

(ii) E[/(a°'-)]^M+(/) ifA>A2. 

Part (i) is a lower bound on the relaxation time and part (ii) is an upper bound. 
Schonmann proved this behavior in dimensions d ^ 2 15j . However, his proof 
left open the question of whether or not Ai = A2. Schonmann and Shlosman 
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settled this question in dimension two, proving that the above holds with Ai = 
A2 ; the transition is sharp in a logarithmic sense |16) . Their proof refines the 
heuristic argument and shows that the critical value of A is a simple function 
of the surface tension of the Wulff shape. The proof takes advantage of specific 
features of the two dimensional Ising model such as duality. When considering 
disordered models, in two and higher dimensions, these simplifying features no 
longer exist. New arguments from the Li-theory of phase coexistence have to 
be used instead. 

The focus of this paper will be the dilute Ising model. For the purpose of 
comparison, we note that the proof of our main result (Theorem 11.2.31 below) 
implies that the upper bound of [16] extends to higher dimensions. We believe 
that our method of proof is valid for all /3 > /3c but we did not make verifying 
this a priority. To avoid certain technicalities we assume that /3 > (3o and (3 ^ J\f 
(see Section [TT^ . 

Let n'^ denote the equilibrium, undiluted Ising measure with a magnetic field 
h > 0. With reference to (j4.2.3l) . the cost of creating a critical droplet under 
E^'^/ft,'^"^ and E^'^ = 0(/3). Define a local observable to be a function that 
only depends on the spins in the region [—1/h, 1/h]'^. 

Theorem 1.1.1. Consider the value 

A, = ^ 



i+1 

Let A > A2. For any positive number Co there is a constant C > such that for 
every local observable f and any h > 0, 



E 



(/fcv.-M)) ^ C'||/||ooexp(-Co//i) 



This is an improvement on the upper bound in |15) and corresponds to 
the upper bound predicted by the heuristic of |16) . Proving rigorously the lower 
bound suggested by [16^ in dimensions three and higher requires the development 
of new arguments which we postpone to a future work. 

The dilute Ising model is a variant of the Ising model that is obtained by 
randomizing the Ising model edge coupling strengths. The impact of dilution 
on the relaxation of Glauber dynamics has been studied in [51 [T3]. In the 
Griffiths phase, which corresponds to the sub-critical regime, the disorder is 
proven to lead to a slowdown of the dynamics. In the phase transition regime, 
the metastability has been investigated for the random field Curie- Weiss model 

We will consider the Ising model on Z'^ diluted in the simplest way possible. 
Independently, delete each edge with probability 1 — p. When p is sufficiently 
large, the remaining edges form a supercritical percolation cluster. From this 
point of view, the Ising model is a special case of the dilute Ising model corre- 
sponding to p — 1. It is natural to ask how the relaxation time depends on p. In 
this paper we show that even a small dilution can greatly reduce the relaxation 
time. 

1.2 The dilute Ising model 

Let Z'' represent the hypercubic lattice. The Ising model assigns each site of Z'' 
a spin of ±1. Let S = {±1}^ denote the set of Ising configurations. 
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Let E — {{x,y} : |lx — y||i = 1} denote the set of nearest neighbor edges 
of Z'^. The equiUbrium Ismg measure with local coupling strengths J = ( J(e) : 
e ^ E) and external magnetic field h is defined using the formal Hamiltonian 

-I E ^(eMxMy) - i ^ M^;), 'TeS. (1.2.1) 

e={x,y}£E xeZ-^ 

We will consider local coupling strengths with the Bernoulli distribution. Let 
Q denote the product measure such that for each edge e, Q( J(e) = 1) = p and 
Q(J(e) = 0) = l-p. 

It is well known that when /i ^ 0, the Ising measure fj,"'-'^ is well defined 
by the Gibbs formalism for any inverse-temperature /3 > and local coupling 
strengths J ^ 0. Consider the spontaneous magnetization of the Ising measure, 

m* = ^lim_Q [m^''' (o-(O))] . (1.2.2) 

When m* > there is said to be phase coexistence. For such f3 there are two 
different Gibbs measures aX h ~ Q, corresponding to the limits h — > 0+ and 
h 0-. 

It is shown in [5] that if the J-positive edges percolate then there is phase 
coexistence in the dilute Ising model at low temperatures. In our settings, this 
means that the critical inverse-temperature 

/3c = inf {/3 > : TO* > 0} 

is finite if and only ii p > pc, where pc is the threshold for bond percolation on 
(Z^i?). 

As well as defining the equilibrium Ising model, the formal Hamiltonian 
defines a dynamic model. Let {a^'~)t ^ o denote a Markov chain on the set E 
of Ising configurations, starting at time with minus spins everywhere, and 
evolving with time according to Glauber dynamics. Given a set of coupling 
strengths J, let Ej denote expectation with respect to the Glauber dynamics. 
Our results extend to some other dynamics such as the Metropolis dynamics 
(see Section [^3)1 . 

A quantity denoted Cdii is defined in SectionlHlthat satisfies Cdii — 0(log jz^) 
as p ^ 1. For the rest of the paper consider p to be fixed in the range (pc, 1). 

Let /?o denote the minimum value such that for all (3 > /3q the assumptions 
of slab percolation (see Section 12. 4p and spatial mixing (see Section 15. 4p hold. 
Let J\f C (0, oo) denote the set of zero measure defined by (|2.2.2p . For the rest 
of the paper the inverse-temperature /3 should be assumed to be greater than 
/?o and not in M. 

For 6 e (0,7r) let denote the cost, up to a factor of h'^^^, of creating a 
critical plus droplet in a cone with angle 9. ~ 0{(39'^~^) and is defined in 
Section [42] 

Theorem 1.2.3. For 9 G (0,7r) consider the value 

^e _ + Cdn9^^ 
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(i) Let A > A^. For any positive number Cq, there are constants C, c > such 
that for any h > 0, for every local observable f , 

^C||/|Uexp(-Co//i)" 
^ 1 - Ccxp(-c/V/I). 

(ii) However close p is to one, at low temperatures the diluted Ising model 
relaxes much more quickly than the corresponding undiluted Ising model; 
with reference to Theorem ] 1 . 1 . 11 as /? — >■ oo, 

-5- inf A^^O. 

A2 ee{o,TT) 

We expect, based on the undiluted Ising model [4], that the slab percolation 
threshold is equal to /3c- Further study of slab percolation and spatial-mixing 
properties for the dilute Ising model would likely extend the domain of validity 
of Theorem 11.2.31 down to the critical point. 

1.3 Heuristic 

The metastability phenomenon for the undiluted Ising model [TB] is related to 
the rate of nucleation of plus droplets with linear size order 1/h. Consider a 
small neighborhood of the origin. Initially all the spins are minuses. Small 
clusters of plus spins quickly form and then disappear. After a short time the 
system looks like it has reached equilibrium with minus spins in the majority. 
However, if we look at a much larger region we will be in for a surprise. A small 
number of larger droplets of plus spin will have formed and started to spread. 
They will eventually merge and cover the whole region, leaving the majority of 
spins in the plus state. 

The rate at which droplets of plus phase form, and what happens to the 
droplets once they have formed, depends on their energy. Let V C M'' with unit 
volume. For 6 > let E^(&) denote, up to a factor of h'^~^, the energy of a plus 
droplet with the shape {b/h)V. E^(&) can be estimated as a balance between 
the surface tension at the phase boundary and the effect of the magnetic field 
h, 

E^{b)/h'^-^ = {b/hf-^T{V) - hm.*{b/hf. (1.3.1) 

Here -^(V) is the surface tension of V (see Section 13. 3p and m* is the mean 
magnetization in the plus phase (11.2.21) . 

Let = The energy function E^(&) is increasing on the interval 

(0, B^) and decreasing beyond . Droplet with b < are unstable and 
tend to be eroded by the surrounding minus spins. Droplets with b > B^ are 
expected to spread. The nucleation of a droplet with b > B^ requires that 
the system overcomes an energy barrier /h'^~^, where E^ := E^{B^). Given 
the inverse-temperature /3 there is a unique shape W known as the Wulff shape 
with minimal surface tension; see the definition of 'yV2n in Section 14.11 Setting 
V — yV minimizes E^. The critical droplet shape is {B^/h)W. 

In any small neighborhood the rate at which copies of the critical droplet 
form is approximately exp(— E^//i''~^). Droplets larger than the critical droplet 
spread out with roughly uniform speed and eventually invade the whole space. 



E 
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The space-time cone of points from which one can reach the origin by time t 
(when growing at a fixed speed) has size 0{t'^'^^). li t = cxp{X/h'^^^) with 

fW 

A > A. ^ (1.3.2) 

then we should expect to see a critical droplet form, and then spread to cover 
the origin, by time t. This heuristic picture has been turned into a rigorous 
proof for the two dimensional Ising model [T^ . 

The dilute Ising model is self averaging so the quenched magnetization and 
the quenched surface tension can unambiguously be defined almost surely with 
respect to the dilution measure Q. It is tempting to try to adapt the previous 
heuristic to the case of the dilute Ising model using the quenched surface tension 
in (|1.3.ip to describe the typical cost of phase coexistence. However, we must be 
careful. In much simpler models, such as random walk in random environment, 
it is well known that a small amount of randomness can change the asymptotic 
behavior. 

Dilution seems to be capable of slowing down the dynamics. Consider an 
expanding droplet of plus phase. If it encounters an area of high dilution it may 
get blocked and have to seep around the obstruction, slowing down its progress. 

However, dilution can also speed up the dynamics. The limiting factor in 
the undiluted Ising model is the rate at which plus droplets nucleate. Nucle- 
ation of plus droplets is infrequent due to the high cost of phase coexistence on 
their boundaries. The dilution creates atypical regions, which we will call cata- 
lysts, where the surface tension is unusually low and so the rate of nucleation is 
unusually high. 

The natural human response to catalysts is to try and classify them. Some 
catalyst do not seem to have much effect on the relaxation time. Consider (when 
d = 2) a circle where all the edges crossing its perimeter have been diluted. If 
a plus droplet forms inside the circle, there is no way for it to spread outwards. 

We therefore want to focus on catalysts that create a sheltered region to 
help plus droplets nucleate, but are not so closed off they prevent plus droplets 
from escaping. There seem to be two competing factors. Large catalysts will be 
relatively rare and so the droplets they help to nucleate will take a long time to 
reach the origin. Conversely, small catalysts cannot do a great deal to increase 
the rate at which critical droplets nucleate. 

We conjecture that there is an optimal catalyst shape that determines the 
relaxation time of the system. However, we do not know how to calculate the 
optimal shape. In this paper we look at a restricted class of catalysts: surfaces 
of diluted edges that form open-bottom cones. We control the nucleation rate in 
the cones, and the subsequent growth of the droplet to regions of more typical 
dilution. This approach leads to an upper bound on the relaxation time that 
is much smaller than the time predicted by the formula (II. 3. 21) with quenched 
surface tension. Indeed, part (ii) of Theorem 11.2.31 shows that asymptotically in 
/3 the values of A2 differ greatly. 

We do not address the issue of the lower bound for the metastable time for 
disordered models. We believe that the more important point is to show the 
existence of the catalyst effect of the disorder. 
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1.4 Outline of the paper 

In Section[2]we define the dilute Ising model and recall some of its basic features. 
The random-cluster representation is used to state a coarse graining property. 

In Sections [3] and |4] we look at the Ising model without a magnetic field. 
In Section [3] we describe the Li -theory of phase coexistence. The theory can 
describe both the typical cost of phase coexistence and the cost of phase coex- 
istence in the neighborhood of catalysts. To combine the two cases we consider 
the cone Ae := {x e M"* : xi ^ ||x||2 cos(6'/2)} where either 6 e (0,7r) or 
6 — 2tt. In Section |4] we look at generalizations of the Wulff shape to Ae- The 
Wulff shape is the shape with minimal surface tension given its volume. The 
Wulff shape can be used to quantify the large deviations of the equilibrium Ising 
model. 

In Section [5] we reintroduce the magnetic field. We justify the energy func- 
tion featured in the heuristic. We prove regularity results concerning cluster 
boundaries. The motivation for this is to study the spectral gap of the dilute 
Ising model in finite regions with various boundary conditions. 

Finally in Section [HI we use the accumulated results to prove Theorem ll.2.31 
We do this by proving that the cone shaped regions act as catalysts. To show 
that the clusters of plus phase formed in the catalysts grow we consider another 
type of cone: space-time cones that are Wulff shaped spatially and growing in 
size with time. 

1.5 Notation 

Throughout the paper C, c, Cstb, Chs, etc, will be used to refer to positive numbers 
that may depend on p, (3 and 9 but not on h. We will recycle C and c to refer to 
various less important positive constants; the values they represent will change 
from appearance to appearance. 

Let S'^~^ denote the set of unit vectors in R''. Let ei, . . . ,6^ denote the 
canonical basis vectors. We will use bold to differentiate continuous variables 
X G R'' from lattice points x lElA. 

Consider ^, S C M'', x e M'' and c ^ 0. Let -h S = {a b : b € yl, b e S} 
denote the sum of the two sets. Let .4 -I- x denote the translation of A by x. 
Let cA denote {ca : a e A}, the set A scaled by a factor of c. 

2 Properties of the dilute Ising model 

2.1 Definition of ^i^^^'^^ 

Let J = (^(e) : e e be a given realization of the coupling strengths. We will 
now define formally the Ising measure /i^'''''' with a magnetic field /i e R and 
boundary conditions € S on a finite domain A C Z'' at inverse-temperature 
;3 > 0. 

Define the external vertex boundary dA of A: 

dA = a+A U d^A where 
d^A ^ {x <^ A:3y e A, {x,y} e E and = ±1}. 
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Taking w to stand for wired, define edge sets for A: 

E{A) = {{x,y}eE:x,yeA}, 
E^{A) = {{x,y} e E : X e A and y e d^A} , 
E'^iA) = E{A)UE^{A). 

The set of spin configurations compatible with C outside A is 

:= {cr e S : Vx ^ A, cr{x) = ((a^)}- 

Changing a Hamiltonian by an additive constant does not change the resulting 
measure; with reference to (jl.2.ip the Ising Hamiltonian H^'^'^ : — )• R can 
be defined by 

^^A^'''(0") = X! •^(e)l{CT(a;)#<7(j;)} + X! ^-^{<^(^) = -l}- 

e={x,y}£E^{K) kSA 

The dilute Ising measure fij^^'^ at inverse-temperature /3 is defined by 

where Z'^''''^ is a normalizing constant, the partition function, defined by 

Zi:^''^= ^ e^p[~(3H';^'^''\a)). (2.1.1) 

We have used a above to index summations over E^. It has also been used as 
a random variable — the mean spin at the origin is written /i^'^'''(cr(0)). Fur- 
thermore, given a set F C Z'' wc will write (t{V) to denote the average spin in 

'^(^) = ME^(^)e (2-1-2) 



2.2 The random-cluster representation for /i^''''* 

The spin-spin correlations in the Ising model can be described by the q — 2 
case of the random-cluster model [TO. We have to be extra careful because of 
the general boundary conditions C S E, dilute coupling strengths (J(e)), and 
the magnetic field h ^ 0. In this section we will describe a random-cluster 
representation for the Ising model n'^''''^ and a joint measure ip'j^'''^. 

The Ising measure ^^^''^ was defined using the graph (A,i?*(A)). Add to 
this graph a ghost vertex g through which the magnetic field will act, and a set 
of ghost edges E^{A) = {{q,x},x G A}. 

When defining the random-cluster model on a given graph, for each edge 
e there is an interaction-strength parameter pe € [0,1]. There is also another 
parameter, q > 0, that influences the number of clusters that are formed. In 
order to describe the correlations of the dilute Ising model we will fix 



2 and 



exp(-/3J(e)), eei?-(A), 
exp{-/3h), eeEs{A). 
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The state space of is f^A = {0, l}S''(A)ui5«(A) ^^1-^ a; G f^A, an edge e is 

open if uj(e) — 1 and closed if w(e) = 0. Two vertices of A are connected if they 
are joined by paths of open edges either 

(i) to each other, or 

(ii) both to a+AU {g}, or 

(iii) both to d~A. 

A cluster is a maximal collection of connected vertices. Let V+, V- denote the 
clusters connected to 9+A U {g}, 9"A, respectively. Let n = n{A,uj) count the 
number of other clusters in A. Label these clusters Vi, . . . , V"„. 

We will define the random-cluster probability measure 4''j^^'^ using a coupling 
probability measure (f'j^''''^ defined on T.'xxflA- The marginal distribution of 
(p'j('''^ on will be the Ising measure iJ.j^^'^. The marginal distribution of 
(fi'j^'^''^ on r^A defines the random-cluster measure (/)a 

With reference to [TDl Section 1.4] define the coupled probability measure 
^p'^pf"^ as follows. Let (cr, S x JIa- Recall the notation (|2.1.2I) for average 
spins. Note that cr{V) = ±1 if and only if a{x) = a{y) for all x,y €V. We will 
say that cr is an w-admissible configuration if 

a{V+)^+l, (t{V-) = -1 and a{V^) = ±1, i = I, . . . ,n. 

If a is oj- admissible, with reference to (j2.1.ip let 



Z 



A 



otherwise let tp'^''''^{{{a,uj)}) — 0. 

For configurations w S Ha under which 9+A is connected to d~A, there are 
no w-admissible configurations. Let C ^Ia represent the set of configurations 
such that the vertices in 9+ A are not connected to the vertices of d~A; is 
the support of 0^'^''*. For lu G D^, there are 2"(^'") oj -admissible configurations 
and 

9n(A.a;) 

' e 

The coupling ip'A^'^ has a probabilistic interpretation. To sample an Ising config- 
uration cr ~ /^a''''* given a sample uj ~ 4'a'''''^, set a{V+) = +1, set (j{V-) = — 1, 
and independently for i = 1, . . . , n set 



1 with probability 1/2, 
- 1 otherwise. 



It is sometimes easier to ignore the ghost edges. Let r(cj) denote the edge 
configuration obtained by closing all the ghost edges. We will say that C A 
is a real cluster if V is an r(a;)-cluster. If F is a real cluster under 0^'''''', but 
not connected to d^A, then 

^^^^ ^1+1 with probability e^''l^l/(l + e^'^l^l), 2 1) 

I —1 otherwise. 
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Let X ^ y denote the event that x and y are in the same real cluster, and let 
A ^ B denote the event that a O 6 for some a G A and b G B. Note that when 
h — 0, all the clusters are real clusters. 

There are two special cases of the h = random-cluster model, wired and 
free boundary conditions. Wired boundary conditions refers to either all-plus 
or all-minus boundary conditions. Under wired boundary conditions = f^A- 
The limit ip"^'"" — lim^^j^d (f)'^^''^ is called the wired random-cluster measure on 
Z"*. Free boundary conditions refers to pretending that 9 A = whilst defining 
^J^QM ^ The resulting measure (t)'''[f'^ depends only on (J(e) : e G i?(A)). The 
measure 4>'^'^ obtained by taking the limit of 0^^'" as A — > Z'' is called the free 
random-cluster measure on . 

Let -f-?- (X) denote the event that the origin is in an infinite real-cluster. The 

set 

{/3 : Q [^"'-f (0 4-> oo)] < Q [Ai^'^'lO ^ oo)] } (2.2.2) 
is at most countable (TTj. It is conjectured that M = 0. 

2.3 Stochastic orderings 

Given two measures /^i,/^2 on a set R^, we will write /ii ^st if there is a 
coupling ((Ti,cr2) on R'^ x R"^ such that 

(i) CTi /ii, 

(ii) (72 ^ and 

(iii) (7 ^ cr' with probability one. 

Holley's inequality [101 Theorem 2.1] can be used to prove stochastic order- 
ings for the ferromagnetic Ising model. The Ising model on a fixed graph A is 
stochastically increasing with respect to the magnetic field and the boundary 
conditions: for hi ^ /12 and any Ci ^ 

The effect of expanding the region depends on the boundary conditions. With 
Ac A, 

fJ-A >st fJ-A but fj.^ Sjst A*A ■ 

Under plus boundary conditions, the random-cluster representation increases 
with h G [0,cx)) and J, 

0a''^'''' «;st 4>'t^''"' if /i2 > /^i > and Ve, J2(e) ^ Ji(e) ^ 0. 

Note that sending J(e) to zero or infinity on the boundary allows us to compare 
free and wired boundary conditions. 

2.4 Coarse graining 

Coarse graining is an important technique in the study of percolation and the 
random-cluster model. The open edges of the random-cluster model percolate 
for f3 > Pc- Slab percolation is a stronger property than percolation [T7]. In 
three and higher dimensions, slab percolation refers to percolation in a slab 
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gd-i X In two dimensions it refers to the existence of spanning 

clusters in rectangles with arbitrarily high aspect ratios. The slab-percolation 
threshold is defined 

/3c = inf{/3 > : slab percolation occurs under ^ /3c. 

It is conjectured that /3c = /3c- 
With K a positive integer, let 

Mk = [~K/2,K/2fnZ'^. 

For i e let B/f (i) := + Ki denote a copy of Mk centered at Ki. 

Let A C Z''. To allow unusually high dilution on the edge boundary of A, 
let J ~ Q and let J' e {0, 1}-^ denote a set of coupling strengths that agrees 
with J in E{A). Let uj E n\ denote an edge configuration for A. 

If BA'(j) C A and, looking at the restriction of w to Bif(i), if there is a 
unique real-cluster A C B/f(i) connecting the 2d faces of Bi<-(«), let B]^(i) = A; 
let B^(i) denote the real A-cluster containing B^(i). Otherwise, let B^(i) — 
Ml,{i) = 0. 

Let Ecg > 0. Recall the definition (jl.2.2p of the spontaneous magnetization 

m* . 

Definition 2.4.1. A 6oa;Bx(«) C A is Ecg-good if: 

(i) B]^(i) is connected (by paths of length one) to eac/i B^(j) such that^Kij) C 
A and 1 1 i — J 1 1 1 — 1 • 

(ii) The diameters of the real-clusters o/ A \ (UjB]^(j)) intersecting Mxii) are 
at most K/2. 

(Hi) B]^(i)nBx(i) contains between K'^m*{l — ecg) and K'^m* {l+Ecg) vertices. 

Otherwise Mji{i) is Ecg-bad. 

Note that any box not entirely contained in A is automatically Scg-bad. If a 
box is 1-bad then it is also Ecg-bad for all Ecg € (0, 1). Recall that we have fixed 
/3 > /3o ^ /3c. The supremum below is over J' e {0, 1}^ that agree with J on 
i?(A). Combining [ITl Theorem 2.1 and Proposition 2.2] yields: 

Proposition 2.4.2. There are constants Ccg ~ Ccg(ecg) > and Kq = /^o(£cg) G 
N such that for K ^ Kq and any Bx(ii), • . • ,Bif (i„) C A, with probability 
1 — exp(— Ccg-fCn), 

sup(y9;( '+'°(BK(ji), . ■ . ,Bi<-(i„) are £cg-bad) ^ exp{-CcgKn). 
J' 

Coarse graining gives a crude measure of the cost of phase coexistence. Con- 
sider a path of neighboring boxes Mxiii), . ■ . , Mxiij) in A. If the first box and 
last box are not connected by a path of open edges, then there must be a 1-bad 
box somewhere along the path of boxes. Moreover, there must be a surface of 
at least [K/Ko{l)\''-~^ 1-bad Bxo(i)-boxes separating Mxiii) from Mxiij) [3 
cf. Lemma 4.2]. We obtain the following corollary to Proposition 12.4.2] Let c^g 
denote a positive constant, independent of Ecg. 
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Corollary 2.4.3. Let K ^ Ka{l). Fork = l,...,n, let Bxiii), ■ ■ ■ ,MK(i'jJ 
denote a simple path of neighboring boxes in A with length jk ^ exp(v^). 
Assume the n chains are disjoint. Let A denote the event that for each k, 
MK{i\) is not connected to Bx(i^j^) in U;^^-^Bx(if )• With Q-probability 1 — 
exp(-4gX'^-in), 

sup(^f expi-c^^K'^-^n). 

We can quantify the extent to which a magnetic field and mixed boundary 
conditions affect the coarse graining property. Recall that £'^(A) denotes the 
set of edges connecting A to the external vertex boundary d^A. 

Lemma 2.4.4. Let C, E Y,. For any B^f (ii), . . . ,Bx(«n) C A, with Q-probability 
1 — exp(— Ccg/^n), 

sup(p^ '^^''(Bx(«i), • . • ,Bi4:(i„) are Ecg-bad) 
J' 

s; exp(/3|£;±(A)| + ph\A\ - c^^Kn). 

Proof. Let (cr, w) be an element of the support of 1^9^''^'''; under w no ghost 
edges are open. Let B denote the set of configurations that agree with (cr, uj) as 
far as the vertices and the real edges are concerned, 

The right-hand side is bounded above by exp(/3|£'^(A)| + /3/i|A|). The lemma 
follows by Proposition 12.4.21 □ 

2.5 The graphical construction of the Glauber dynamics 

For ^ G E^, let {o'f'^)t^s denote the dynamic Ising model, started at time 
s in state ^ and evolving according to the Glauber dynamics (also knows as 
heat-bath dynamics). The Glauber dynamics can be described by the following 
graphical construction. Let (t|'^ = ^. Place a rate-one Poisson process at each 
vertex. Label the points of the Poisson process {xi,ti) with ii < t2 < • • • ; to 
each point {xi,ti) attach a uniform [0, 1] random variable Ui. Let (Tfl5 denote 
the Ising configuration immediately before time t. For each point of the Poisson 
process, resample the spin at Xi from the Ising measure conditional on the 
state of the neighboring spins. If the probability cr^Xi) = +1 conditional on 
{cr(?/) = (J(.'l(?/) : y ^ x} is q, set alf{xi) = +1 if ?7i > 1 — g, and —1 otherwise. 
The dynamics are: 

(i) Monotonic, if ^ < 77 then CTj'^ ^ ctj'''. 

(ii) Finite range, to update site x only requires knowledge of the neighbors. 

(iii) Bounded with respect to the transition rates. 

Our results are also valid for other dynamics, such as the Metropolis dynamics, 
that share these properties. 



12 



3 Li-theory 



3.1 Microscopic and mesoscopic scales 

Recall that we have fixed p £ (pcA) and /3 > /3o with (3 ^ N. 

To take advantage of the coarse graining result, we introduce some notation. 
With reference to Theorem [1.2.3l let h> 0. Define 

if = [/i-i/^^-i)], N = K[h-^/K\Kh-^. (3.1.1) 

We will call A'^ the macroscopic scale. This is the scale at which nucleation 
of plus droplets occurs. The number K denotes a mesoscopic scale. We will 
consider regions with size order N composed of boxes Bi<-(i). 

Let I? C M'' denote a connected region. Let D(I?, N, K) denote a discretized 
version of T) composed of mesoscopic boxes, 



K K'"^ 
'2N'2N 



CV}. (3.1.2) 



We have used h > to define a set A = B{V,N,K) with size 0{N) ^ 0{h-^) 
on which we wish to study the Ising model Given A, we will also want 

to consider the Ising measure /^a^'". From now on, h will always determine 
the scale A^ but will not always indicate the strength of the magnetic field. 
The coarse graining implies that under /iA^'", (7(8^^ («)) is close to either ±m* 
with high probability. This motivates the definition below of the magnetization 
profile : R'^ — > M associated with the Ising configuration a. 

Let A denote an arbitrary finite subset of Z'^; we will mostly be interested 
in sets of the form 'U{'D, N, K) but we will also consider sets that differ from 
D(2?, A^, K) around the boundary. Let C G S denote a boundary condition. 
Recall the notation (j2.1.2p and that under ^I'jy'"^ , <t{x) = ({x) for x ^ A. 

For cr e E^, define the profile : R'^ ^ R as follows. For x e R'^, choose 
i such that A^x e [-K/2, K/2Y + Ki. Let 




c7(IBA-(j))/m*, IBA-(i)cA, 
o-(BA'(i)), otherwise. 



(3.1.3) 



A value of 1 indicates plus phase, indicates minus phase. The idea behind 
Li-theory is that can be approximated by the class of bounded variation 
profiles. The large deviations of can be described in terms of surface tension. 



3.2 Surface tension in a parallelepiped 

In statistical physics, surface tension is the excess free energy per unit area due 
to the presence of an interface. The definitions of surface tension in [TO] and 
[l8] differ by a factor of (3; we have chosen to follow [18]. 

Let (n, U2, . . . , Ud) denote an orthonormal basis for R'' and let TZ denote the 
rectangular parallelepiped 

'RL,H{-n., U2, . . . , vid) < tin + ^ tkVik ■ t € 

I k=2 



H 


H 


X 


L 


L 
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TZ is centered at the origin, has height H in the direction n and extension L in 
the other directions. 

Let A = D(7?., iV^ 1) denote a discrete version of A p.l.2p . The box A has 
sides of length NL in the directions U2, . . . , u^. The surface tension can be writ- 
ten in terms of either the Ising model partition function (j2.1.ip or the random- 
cluster representation. 

Definition 3.2.1. Let C, denote the configuration in S given by C,{y) — if 
y ■ n ^ and C,{y) = —1 otherwise. The surface tension is defined by 



{NL)d-i ^^^.C,o (^NLY-i %J^+'0(^Biy 

Let J ^ Q. Surface tension converges in probability as — > 00 [TH Theo- 
rem 1.3]: 

Proposition 3.2.2. For /3 > and n S S'^~^, there exists T(n) ^ 0, the surface 
tension perpendicular to n, such that for all parallelepipeds TZ — TZl,h{'o., U2, • . • , u^), 

Tj{ !• T(n) as N ^ 00. (3.2.3) 

Surface tension is strictly positive at temperatures below the threshold for 
slab percolation [TBI Proposition 2.11]: 

Proposition 3.2.4. There are constants C, c > such that for n E S'^^^ and 
13 > /3c, cT(ei) ^ T(n) < CT(ei) ^ C/S. 

We note for completeness that we are discussing quenched surface tension. 
Annealed surface tension, which will not be used in this paper, describes the cost 
of phase coexistence under the averaged measure Q/i''^'^'''. When studying large 
deviations under the annealed measure, the environment J changes to reduce 
the surface tension. Although we are interested in the large deviations of J, we 
prefer to control them 'by hand' using the Qg notation defined in (j4.3.ip . This 
is less efficient in terms of the size of the large deviation needed. However, it is 
much simpler. 



3.3 Surface tension in cones 

Li-theory describes the Ising model at equilibrium [51171118). We will restrict our 
attention to certain subsets of R'' with zero surface tension at the boundary. At 
the microscopic scale, this corresponds to the sampling J ^ Q but conditioned 
on the existence of a surface of edges with J(e) = 0. 
For 9 G (0,27r] define a linear cone Ae, 

Ag :- {xeR'^-.xi^ ||x||2 cos(0/2)}. 

For 9 — 27r, Ag is simply the whole of M.'^. Let denote the Lebesgue measure 
on Ae and let V{u,e) denote the £''-ball of radius e about u : Ag ^ M. 

The perimeter V{U) of a Borel subset U C Ag can be written in terms of 
functions of bounded variation (see [2 Chapter 3] and (TS] Section 3.1]). The 
set of bounded variation profiles BV is given by 

BY -.^{U -.U C Ag is a Borel set and P(C/) < 00} . 
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Bounded variation profiles U G BV have a reduced boundary d*U and an outer 
normal n'^ : d*U — > 5"^^^. Let "W*"^ denote the d — 1 dimensional HausdorfF 
measure on Ae\ 'H'^^^{d*U) — 'P{U). We will write dU to refer to the reduced 
boundary of U excluding (when 9 < 2tt) the boundary of Ae , 

dU ■.= d*U\d*Ae. (3.3.1) 

The outer normal n'^ defined on dU is Borel measurable. With reference to 
(13.2.31) , this allows us to define the surface tension and energy of bounded vari- 
ation profiles for the dilute Ising model. Define the surface tension by 

F{U)^ I r(n^(x))d-H''-i(x), [/ e BV. (3.3.2) 

JdU 

Define the energy E by 

£{U)^F{U)~ l3m*C'^{U), C/ G BV. (3.3.3) 

The motivation for these quantities is that they measure, in the following sense, 
the cost of phase coexistence associated with the Ising model. Sample J from 
Q conditional on J{e) — for all e = {x, such that x but not y is in 
^{Ag,N,K) p.l.2p . This makes the surface tension on d*Ae zero. Let T) 
denote a compact subset of Ae- For a profile of bounded variation U C V, the 
surface o{D{U,N,K) has size Oil/h"^-^) and D{U,N,K) has volume 0{l/h'^) 
(|3.1.ip . Heuristically, we expect that the probability of seeing the plus phase in 
D(C/, N, K) and the minus phase in D(I? \ U, N, K) to be approximately 

J,-,h 

The infimum is over profiles U' € BV compatible with the boundary conditions. 
For general D, it is difficult to evaluate the infimum — the conflicting contri- 
butions of the positive field and negative boundary conditions may lead to a 
complicated equilibrium magnetization profile under Mjj'fi^'Ar k)' "^^^ problem 
is simpler if T) is the Wulff shape. 

4 The WulfF shape in A9 

4.1 Wulff, Winterbottom and Summertop shapes 

Let U C Ae denote a set of bounded variation. Consider the problem of mini- 
mizing T{U) given that U has volume 

Proposition 4.1.1. Let 9 £ (0,7r] U {27r}. The problem of finding a set of 
bounded variation U C Ae with volume b'^ and minimal surface tension has a 
unique solution when 9 < n; for 9 — n and 9 = 2tt the solution is unique up to 
translations. There is a scaling constant we such that the solution is the convex 
shape 

We{b) = web{^ e A : Vn e S'^-\ x • n < T(n)}. (4.1.2) 



exp 



HU) \ 



exp 



■mf£iU')-£iU) 



under 
under 
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If we omit the b, take b = Wg^ so that We = Weiwg^). 

Special cases of We are known by a variety of names. When 9 — 27r, We is 
the Wulff shape. When 9 — tt, We is the Winterbottom shape. When 9 S (0, tt) 
and d — 2, We is the Summertop shape [H]. We will refer to We as the Wulff 
shape in Ae- 

Proof of Proposition \4.1-1\ Let U denote a compact subset of Ae ■ With refer- 
ence to 11, (G)], the formula p.3.2p that defines the surface tension J-{U) is 
equivalent to 

^ \iU + eW,.)C^Ae\-\U\^ 

e^O £ 

We can use (j4.1.2p to define a second measure of surface tension. Let 

PiU) ^ hm l^ + ^^-'-'^l 

e^O e 

Observe that: 

(i) For any x € y^e and £ > 0, 

X + eWe C (x + eW2-K) n Ae- 

Thus for U e BV, P{U) < T{U). 

ill) By the Brunn-Minkowski theorem, We(b) is the unique shape in Ae (up 
to translations) with volume b'^ and minimal J^-surface tension. 

(iii) By the convexity of W2-K, when U = We, 

{U + eW2n)f^Ae = U + eWe 
and so P{We{b)) = F{We{b)). 

Therefore any shape with volume b'^ in Ae with minimal J-'-surface tension must 
take the shape We{b). The claim of uniqueness when 9 < tt follows from the 
fact that for any x e A \ {0}, J"(x + We{b)) > T{Weib)). □ 



4.2 Critical droplets 

Let E^(fe) account for the cost of filling We{b) with the plus phase, 

E^{b) := £iWe{b)) = b'^-^T{We{l)) - b'^^m*. (4.2.1) 

The positive term represents the cost of phase coexistence, the negative term 
represents the benefit of conforming to the magnetic field. Let denote the 
maximizer of E^, and let Sfoot denote the positive root of E^, 

, d^l ^iWejl)) e nmil)) 

The significance of i?foot that the Ising measure with minus boundary condi- 
tions and magnetic field h on I]>{We{b), N, K) favors the minus phase if < -Bfooti 
whereas it favors the plus phase if 6 > B^^^^ . Let 

(T{We{l))Y fd-1^ ''-^ 



EW)=(^^^^^) (^) . (4.2.3) 



16 



The maximum of E characterizes the energy needed to create arbitrarily 
large plus droplets in the cone Ae, starting from the minus phase. 

As the J(e) are independent we can find regions that resemble, due to high 
local dilution, Bi(We{b), N, K) for any 9 £ (0, tt) and any & > 0. However, to 
maximize the number of catalysts we do not want to take b any larger than we 
have to. 

Proposition 4.2.4. The diameter of the critical droplet is bounded uniformly 
over I3> ande e (0,7r) U {27r}. As 9 ^ 0, E^ = 0(139'^-^). 

Proof For U g BV, J"(f/) has order (SW^-^idU) by Proposition Choose 
ug such that the cone 

1^9 we{x e Ae : xi 1} 

has unit volume. As — > 0, both Ug and Wg{l) have length of order 0-i<^-^)/<^ 
in the xi -direction. By the optimality of the Wulff shape, FiyVe{^)) =^ J^iUe) 
which has order (3{9ue)'^~^. Substitute this approximation into (|4.2.2p and 

mM- □ 
4.3 Notation for Wulff shapes 

Consider the discrete analogue Ag of Ag at microscopic scale N and mesoscopic 
scale K ((3T2)) . 

Ag ■.= B{Ae,N,K), 
and the discrete analogues of the Wulff shape, 

We(6) :=]D)(We(6),iV,i4:), b^O. 

For 9 ^ 2tt, the edge boundary of Ae is infinite, thus the probability of finding 
a pattern of dilution that carves out a translation of the Ag anywhere in Z'' 
is zero. Instead let B^ax > denote a fixed, but as yet unknown, quantity. 
We will limit our attention to the region We{B^^^). To impose free boundary 
conditions on the portion of the boundary of Wg (i?max) corresponding to d*Ae , 
let Qe denote the dilution measure Q conditioned appropriately, 

Qe := Q[ • \ yx^y such that x e Wg{B^„,,J but y ^ Ag, J{{x,y}) = 0]. 

(4.3.1) 

The probability of seeing such a pattern of dilution is simply (1 — p) raised 
to the power of the number of edges that are conditioned to be closed in the 
definition of <Qg. Wg{b) has size order along the xi-direction and 

order b9^^'^ in the directions X2, - ■ ■ ,Xd- The number of edges that need to be 
diluted is therefore order (B^^J"^'^-^''/'^) x {B^^J^^'^)'^-'^ ■ The Q-probability 
of the event conditioned on in (|4.3.ip is thus 

exp (- {bU)'" 0-"'0 (log ^) /h'^-') . (4.3.2) 

As well as Wulff shaped regions, we also need to consider Wulff 'annuli': the 
difference between two Wulff shapes. With ^ 6i ^ 62 =^ ^max: I'^f 

Wg{biM)-^m{b2)\Wg{bi) and W(,(6i, ^2) We(62) \ We(6i). 

When 9 = 27r, yVg{bi,b2) is annular; otherwise it is simply-connected. There 
can be three parts to the boundary of 'Wg{bi, 62): 
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Figure 1: From left to right: We is the intersection of Ae with W2ir- The set 
We(5i,62) with (+,—), (— ,+) and (+,+) boundary conditions. The dotted 
hues indicate free boundary conditions. 

(i) the inner boundary dWe{bi), 

(ii) the outer boundary 9 We (^2), and 

(iii) the free part of the boundary d*Ae n 9*yVe(&i, ^2). 

If &i = then there is no inner boundary and W6/(6i,&2) = ^eib-z)- It = 2tt 
then the free part of the boundary is empty. 

Let (+, — ) denote an Ising configuration that is equal to +1 on We{bi), and 
equal to —1 on Ag\ We (62); see Figure [TJ We will show in Section that 
for b € [&i,&2], the probability that M^'"-* is close to yVe{b) under is 
approximately 

f :F{W9ibi))-:F{We{b)) \ 

h^^ )■ 

As well as (+, — ) boundary conditions, we will also consider boundary conditions 
of (— , +), (+, +) and (—,—). We may simplify (+, +) to + and (— , — ) to — . 

We will also need to consider some sets that only differ from Wg{b) and 
We(6i,&2) at the mesoscopic scale. The sets Wg{b) are subsets of the discrete 
cone Ag. Let {xi, X2, ■ ■ ■} denote an ordering of Ag and let 

-{xi,...,x„}, n^O. (4.3.3) 

As a function of 6, |W6((&)| is non-decreasing. 'Wg{b) is composed of boxes of 
K'^ vertices, so \'Wg{ • )| is a step function (i.e. piece- wise constant and cadlag). 
We can assume that the ordering has been chosen so that for 6^0, Wg{b) ~ 
A^J^o'-''^^ Let W'g{-) denote a second family of increasing subsets of Ag such 
that 

(i) Wg{b) = Wg{b) at the points b of discontinuity of \Wg{b)\, 

(ii) Wg{b) C W'sib), 

(iii) for each n, for some bg{n), Wg{bg{n)) — Ag. 
Let W^(6i,fe2) - W',(fe2)\We(fei). 

4.4 High and uniformly high Qe-probabihty 

We will say that an event occurs with high Qg -probability if under Qg it occurs 
with probability at least 1 — C exp{—c/ Vh). Note that by taking C large, we 
only have to consider h £ (0, ho) where ho can be arbitrarily small. 
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Given B^^^^ > 0, suppose that 

61 s$ 6 ^ 62 ^ S^,ax and A = W;,(6i,62). (4.4.1) 

Given a class of events defined in terms of b,bi and 62 we will say that they 
occur with uniformly high Qg -probability if each event occurs with probability 
at least 1 — C exp{—c/ Vh), uniformly over (|4.4.ip . Abusing this notation, we 
may place some additional restriction on bi and 62; for example fixing bi — 0. 
In that case interpret (|4.4.1[) with the additional restriction in place. 



4.5 Li-theory under (w, — ) boundary conditions 

The Li-theory developed in [18. describes the dilute Ising model in cubes A = 
{1,...,A^}'' under the measure J ~ Q. The proofs in [18] are easily 

adapted to sets of the form Wg(6i, 62) with J ^ Qg. Moreover, the methodology 
accommodates the (w, — ) boundary conditions described below. 

Consider A as in (|4.4.ip . Let (w,— ) denote wired boundary conditions on 
the inner boundary of A, and minus boundary conditions on the outer boundary 
of A. If bi — 0, (w, — ) simply means minus boundary conditions. 

This is equivalent to starting from (+,— ) or (— ,— ) boundary conditions, 
and then replacing the inner boundary of A with a single Ising spin variable. 



J,(w, — ),0/ r 1 \ 




This is a measure on S^^' •* U ' •* . Let 9" A denote the wired inner-boundary 
of A. Let ct(9'"A) = ±1 according to whether a e ^^{^'^^ or e S*^"'^'. Define 
to be either or M^^'"^ according to (j{&"k). 

The (w, — ) measure with /i = has a natural random-cluster representa- 
tion with wired-inner and wired-outer boundary conditions. The corresponding 
coupled measure is 



¥>A TjmTo ■ (4.5.1) 



Note that only the s = a-{d^A) term in the numerator of (|4.5.1I) is positive. 

Let Dj^'"-* refer to the event that the inner boundary d^A and the outer 
boundary d~A are not connected in the random-cluster representation. For 
e > 0, let Db g denote the event that the inner-boundary takes the plus spin and 
that the phase profile is in a neighborhood of yVe(6), 

D,,, Z?i"-' n {a(a-A) = 1} n {m^T'-) e V(We(6),£)} . 

Parts (i) and (ii) below follow from the proofs in [TB] of Proposition 3.10 and 
Theorem 1.11 respectively. 

Proposition 4.5.2. Let b G [61, ^2] and Evim > 0. With high Qe -probability: 
(i) The probability of phase coexistence is bounded below, 

/z'^-ilog^:('("^")'°(D,,,_) ^ -T{We{b))-e^^. 
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(ii) The probability of phase coexistence is bounded above, 

4.6 Large deviations under (+, — ) boundary conditions 

With b and A as in (|4.4.ip . we will give upper and lower bounds for the cost 
of phase coexistence under mixed boundary conditions in the absence of an 
external magnetic field. 

Under (+, — ) boundary conditions, the Ising measure favors the minus phase 
because the minus boundary is bigger. Let > 0. Large deviations of the 
magnetization defined in p.l.3p away from the minus phase are con- 

trolled as follows. 

Proposition 4.6.1 (Upper bound). With high Qg -probability 

h^-Hogf^'/^'-^'" (^J M^+^-KlC" ^ b''^ ^ T{We{bi)) - T{Wg{b)) + Spm- 

Proposition 4.6.2 (Lower bound). With high Qg -probability 

h^-Hog//+'-^'" (^J M^p-^ dC' ^ b""^ F{Wg{b^)) - F{We{b)) - e^r^. 

Proof of Proposition \J. 6. 1\ By the definition of the (w, — ) measure, 



^ ( / m(^'-' dC'-^^b''^ -^ ^ 



■/,(+,-),o 



Proposition 14.5.21 part (ii) provides an upper bound on the first term on the 
right-hand side of (|4.6.3p . 

and s = — 1 terms in the numei 

side of (|4.5.ip are equal, so 



For ui € D^^' the s — and s = —1 terms in the numerator of the right-hand 



Summing over w, 



, J,(w -),0/ X _ ^^A 0A {^) 



^.,w,-),0(i3w.-) ^ (4.6.4) 

By Proposition 14.5.21 part (i) with b = bi, 
'/'A (Da )>fA (Dfci.e™) ^ exp ( I. (4.6.5) 
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Combining inequalities (|4.6.4p and (|4.6.5p produces an upper bound for the 
second term on the right-hand side of (|4.6.3p . The proposition fohows by taking 



Proof of Proposition \4.6.S\ With £wm > let b' = {/¥ + Ewm so that Df,/^e„^ 
implies / M^'"-* dC^ ^ fe'' . By equation (|4.6.4p and the definition of the (w, -) 
measure, 

n J,(W. — ),0 /-r^ \ 

Proposition 14.5.2] part (i) gives a lower bound for the numerator on the right- 
hand side. Proposition 14. 5 .21 part (ii) with b = bi gives an upper bound on the 
denominator on the right-hand side. Take e^m sufficiently small. □ 

Proposition 4.6.6. With reference to (|4.4.ip . for fixed Spm both Proposition \4.6. 1\ 
and Proposition \4.6.'^ hold with uniformly high <Qg -probability. 

Proof. Consider Proposition 14.6.11 the case of Proposition 14.6.21 follows sim- 
ilarly. Controlling uniformly the Qg-probability can be reduced to the con- 
trol of a finite number of events. Let AI{bi,b2,b) denote the increasing event 

|cr : / M^'"^ dC^ ^ fe'*! with M^'"^ defined with respect to W^(&i,62)- 

The measure /i^''^' '' '' is stochastically increasing with bi and 62- Let C, e > 
and let 



b[=e\b,/e^, &'2=e[&2/£l, b' ^ e 



e~'{bd-Ce{BUr-^ 



For triples (61, 62, 6) satisfying (|4.4.ip . the triple {b[, b'2,b') takes a finite number 
of values. With high Qe-pi'obability, 

/.'^-iiog^4f|-g {M{b',,b'„b')) ^ Hmm-Hm{b')) + ^. 

C can be chosen such that li F{'We{bi))—J-iyVe{b))+e-pi-a < and e is sufficiently 
small then M(6i, 62, &) =^ Af 62, and 
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5 Spatial and Markov chain mixing 

We will describe the dilute Ising model at equilibrium under mixed boundary 
conditions and a magnetic field. These results will be used in the next section 
to show that sufficiently large plus droplets spread in a predictable way. 



5.1 Stability of minimum energy profiles 

Consider the Ising measure with a magnetic field h and (+,— ) boundary con- 
ditions. Recall the definitions of the energy functions 8 (|3.3.3p and (|4.2.ip . 
With reference to (j4.4.ip . consider the case E^(6i) > E^(52). The minimum 
value of E^(6) is attained when b — b2. Geometrically, this means that the 
profile U with We(&i) C U C ^0(62) that minimizes £{U) is Weib2). The 
minimizer is unique and stable. 
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Proposition 5.1.1. Let £stb G (0, 1) such that bi < 62(1 — £stb) and E^(&i) > 
£^(62(1 — Estb))- There is a constant Csth — Cstb(estb) > 0, independent of B^^^, 
such that for profiles U S BV, 

C%U) e [bf, 6^(1 - e,tb)1 =^ £{U) ^ £{We(b2)) + 262Cstb. (5.1.2) 

Given Ssth, with uniformly high Qg-probability, 

(/ M^^-^ dC' ^ bin - e.,,r'^ ^ exp . (5.1.3) 

Proof. We will first show inequality (|5.1.2p with 



Cstb := . ,9 dr. — — where e 1 - (1 - e^tb) • 

By the optimality of the Wulff shape, Proposition I4.1.1[ we can assume that 
U = Weib) for some b. The unique maximum of occurs at B^. As E^(6i) > 
E^(62(l — Estb)) the function E^ must be decreasing in the region [62(1 — £stb), ^2]; 
it is optimal to consider 

b ^ b2{l - Sstb) > Bl (5.1.4) 

With b as above, 

E^(&) - E«(&2) = TiWeimb"-' - bf) ~ Pm*{b^ - b^) (5.1.5) 
= T{We{l))b't\{l - ef''-^^''' - 1) + I3m*ebi 

By and ([5Xi| . 



Let /(e) = (1 - e) - For e E (0, 1), /'"(e) > so 

/(e) - /(O) - e/'(0) ^ /"(0)eV2, e e [0, 1], 
or more explicitly 

(l_e)_^T^ + £l_i^£2.^^ ee[0,l]. (5.1.7) 
Substituting (j5.1.6p into (|5.1.5D and then using (|5.1.7p gives 

E'{b) - E«(62) > -F(W,(l))6^i . . ^ 262Cstb 

as required. 

We will now show inequality (j5.1.3p . Let S count, up to an additive constant, 
the number of plus spins in A, 



cr{x) + m* 
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By the definition of A, M^t' ' and S, 



-.*hd I l,d(. 



m*bt + h'^S - m* 



0{K/N). 



The magnetic field corresponds to a Radon-Nikodym derivative controlled by 
S. With Z := [exp(/3/iS')], 

We will find a lower bound on Z. When h is sufficiently small, 
Applying Proposition 14.6.21 with b — b2, 

h"-^ logZ ^ /i^-i logMl''+^"^-°(M(+'-) e V(We(&2),epm)) + /3m* - b^) - 0(e 

> ^(We(6i)) - J-(We(62)) + /3m*(fe^ - 6^) - O(ep„0 

> £{yVe{bi)) - £{yVe{b2)) ~ 0{e^^). 

Let n = [(&2(1 - Estb) - 6i)/epmJ- We can write 



u 



be{bi,6i+epm,...,6i+epmn} 

By Proposition 14.6. 11 for b e &i + epm, . . . , 6i + epm'^}, 

/.^(+'-''° exp(/3/i5); ^ M^^t'-) d/:'' e [6^ (6 + £p,„)^] 

^ exp (/i^-'^ - + /3m* [(6 + e^,,-,f - bf] + 0(epnO]) 

< exp {h^-'' [£{We{bi)) - £{We{b)) + 0(£p,„)]) . 

The left-hand side of (|5.1.3p is therefore at most 

exp {h^-''[£{We{b2)) - £{m{b)) + O(£p„0]) • 

be{bi,6i+epm,---,bl+epmn} 

Taking Cpm small with respect to Sgth, inequality (|5.1.3p follows by (|5.1.2p . □ 

Consider now the case E^(6i) < £^(62)- The optimum profile matching 
(+,— ) boundary conditions is yVe(fei) so the minus phase is dominant. 

Proposition 5.1.8. Let £stb > 0. Suppose that {bi + £stb) < £^(62)- There 
is a constant c'^^y^ — Cgj)^(estb) > such that with uniformly high Qg -probability 



JA+-},h 

Ma 



(^J U^+^-^ dC' ^ (61 + e^thf^ < exp (^-^^ ■ 
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We will omit the proof of Proposition 15.1.81 as it is similar to the proof of 
Proposition 15.1.11 

Now consider the case 61 = and 62 > under minus boundary conditions. 
In order to get a stability property that does not depend on the sign of £^(62) 
we will condition on seeing a large region of the plus-phase. With ec > let 

fl''^'^^ := \C) where (5.1.9) 

C = C{ec) := I a e S;, : / d£^ ^ [Bl f I . 

"'We(S«+ec) J 

Proposition 5.1.10. Let estb,ec > 0. Suppose 61 = and 62(1 — Estb) > 
Bf^ + ec ■ Recall Cgtb from Proposition \5.1.1\ Given Ssth o,nd sc , with uniformly 
high Qe -probability, 



exp 



^2Cstb 



Proof. The proof of (|5.1.2p (with B^ playing the role of hi) implies that for all 
U eBY, 

£'*(£/) G [(Bf)^ 6^(1- e.tb)1 =^ £iU)^£{We{b2)) + 2b2c,tb. 

Therefore £"(62(1 - £sth)) - £^62) ^ 2b2Cstb . 

Returning to the context of Proposition 15.1.101 we have bi = 0. Let a = 
min{0, £^(62)} denote the minimum of E^(6) for 6 G [61,62]. Let b denote the 
minimum value in the range [62(1 — £stb), 62] such that 

bi-b'^^ {Bl + ecY - [BlY . 

Treating the magnetic field as a Radon- Nikodym derivative as in the proof of 
Proposition 15.1.11 



fl~E^(b2)-3b2Cstb/2 



and 

'a_E»(62)-62Cstb/2 



□ 



We will now consider two different boundary conditions. By ProDOsition l4.6.11 
plus/minus symmetry when /i = 0, and monotonicity, the plus phase is dominant 

l^A 



under 



Proposition 5.1.11. Let Ssth > 0. There is a constant c"j^ = c'^^^^{esth) > 
such that with uniformly high Qe -probability, 



l^A 



[ m(^-'+) dC ^bi- (61 + £,tb)' ) exp (-^) 



Finally, consider boundary conditions of plus on the inner boundary and free 
on the outer boundary. 
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Proposition 5.1.12. Let £stb > '^^^ suppose 62 ^ 2£stb- With uniformly high 
Qg -probability, 



We omit the proof as it is similar to the others in this section. 
5.2 Phase labels in 

Lemma 12.4.41 provides a simple measure of the cost of phase coexistence in a 
region conditional on the boundary. Using phase labels — defined in terms of 
the coarse graining — to describe the boundary conditions allows for a sharper 
bound. 

Take A as in (|4.4.ip . Let lo e JIa, and let a E denote an oj-admissible 
spin configuration. Define phase labels in terms of the coarse graining with 
£cg = 1: 

-1, Mxii) is 1-good and cr{M^j^{i)) = +1, 
= { -1, ^K{i) is 1-good and cr(Bi(i)) = -1, (5.2.1) 
0, otherwise. 

Let F — Uig7Bif(i) denote a subset of A composed of whole mesoscopic boxes. 
Let dl denote the set of points i e Z'^ at ioo-distance 2 from / such that 
V>K{i) C Ag. Let ij} : I dl ^ {±1,0} denote a label configuration assigning 
labels to the boxes IB/<-(i). Let Int^ (respectively Ext^) denote the event that 
^(i) = ?A(i) for i e / (respectively dl). Given -0, let fj count the number of 
labels in /. Let fgi,fgj,fgj count the number of boxes in dl with phase 
labels +1, —1 and 0. Let count the maximal number of disjoint paths (the 
maximal flow) inside / between the +1 labels and —1 labels. 

We can find constants Ccut , Ccut > such that the following hold. 

Lemma 5.2.2. With Qg -probability 1 - exp(-Ccut/r^''"^)' 

ipi:^'''{lnt^ n Ext^) s; exp {C^tifgi + fgi)K^'' - c^mfT K''-^) ■ 

Lemma 5.2.3. With Q^g -probability 1 - cxp(-Ccut[.//i^ + ff^K'^^^]), 

n Ext^) exp (Ccut min{/a", + /g, , fgj + /+ }if ''-^ 

+ph\T\A'' - c,^t[fiK + frK'-'])- 

Proof of Lemma \5.2.2\ Let Ffe denote the union of the mesoscopic boxes at dis- 
tance at most k from F, 

rfe=|J{B/f(j):3«e/, ||«-j||oo ^ fc}. 

Note that F C Fi C F2, and for i e dL, B/f (j) C F2 \ Fi. Write w = Wg ® uj^^t © 
Wext where is the configuration of the ghost edges and Wint is the configuration 
of the real edges £'(Fi). 

We would like to condition on the event Ext^. However, the resulting mea- 
sure is complicated because Ext^ carries information about not just Wcxt and 
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(cr(B^(i)) : i G a/), but also about Wint- Let Ext^ denote the event that ojcxt, 
and the states of the vertices in A \ Fi, are compatible with Ext^. Let Int^ 
denote the event that Wjnt is compatible with Int^. 

The coupled measure (f'j^'''''' has a property related to the finite energy prop- 
erty of the regular random-cluster model. The probability of edge e being closed, 
conditional on all the other edge states and all the spin states, is bounded away 
from 0: 

ip-j:^'''He) = I a and (w(/)) f^,) ^ 1 - p,. (5.2.4) 

This tells us the cost of conditioning on edges being closed. Surgically closing 
edges saves us from having to consider mixed boundary conditions. 

Let T denote the event that all edges spanning between a box Bx(j) with 
i € dl and ^{i) ^ and a box B^O) C Fi are closed. By (|5.2.4I) . for some 
constant Ccut > 0, 

'"(Int; n Ext:^) '"(Int:^ n Ext; n T) exp (Ccut(/a"/ + fai)K''-^) ■ 

(5.2.5) 

Suppose that cjcxt splits A \ Fi into ojext-clusters W^+, Wi, . . . , Wn- We 
stress that Wext is a partial edge configuration: the Wext-clusters in A \ Fi may 
be connected by an (wg © a;int)-path. Assume that the event Ext^ n T holds. If 

Ml^{i), i e dl, intersects Fi then ^{i) = 1. Without loss of generality we can 
assume that 

(i) Wi, . . . ,Wi correspond to B]^,(i) with i e dl and — +1, 

(ii) VFi+i, . . . , Wm correspond to clusters with diameter less than K/2, 
(in) W,n+i, . . . , Wn (and W-) are not connected to Fi. 

Consider the conditional measure ipj^'''^{- \ Ext^ HT, Woxt)- The clusters W+ 
and Wi,...,Wi act as plus boundary conditions. The spins of the clusters 
VFi+i, . . . , Wm are unknown so the clusters act as wired boundary conditions. 
Let c^'* denote the marginal measure on cjjnt corresponding to Lpj^'^''^{ ■ \ Ext^ n 

^^'^■''(Int^nExt^nr) ^ sup 4>Hlnt'^)- (5.2.6) 

WextSExt^nT 

Let A denote the decreasing event that there are no Wint-open paths in Fi 
between B^^-boxes with = +1 and — —1; note that Int^ C A. By 
monotonicity, as A is a decreasing event and cj)^ is increasing with h E [0,cxi), 

4>\lnt'^) ^ 4>'\A) ^ 00(A). (5.2.7) 

With reference to Corollarv l2.4.3l we can find a collection of //^ disjoint chains 
of boxes such that A implies that for each chain, the first box is not connected 
to the last box. With Qg-probability 1 - ex.p{-c'^^ff^ K'^''^), 

^°{A) < eM-<jrK'-'). (5.2.8) 
Collecting (|5.2.5p - (|5.2.8p gives the lemma. □ 
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Proof of Lemma \5.2.3\ We will first consider the case /g} ^ f^j. Let Fi, 
Int^, Ext^, T, Wi,. . . , Wn and be defined as above. 

Recall inequality (|5.2.6p . With reference to (j2.2.ip . the sizes of the Wi+i, . . . , W„, 
affect the interaction (under 0'') of open clusters in Fi with the external mag- 
netic field. Recall that the clusters Wi+i, . . . , Wm have diameter at most K/2: 

4>''ilnt'^)/4>%lrit'^) < exp(/3/i|Fi U Wi+i U • • • U s=: exp{/3h\T\A''). (5.2.9) 

By Proposition l2.4.2l and Corollarv l2.4.3l there is a positive constant c > such 
that with Qe-probability 1 - exp(-c max{/}'X, ff^K''-^}), 

0"(lnt'^) sc: exp(-c max{/OA', frK^-^})- (5.2.10) 

Collecting (|H:^ - (IHXH11 and (I^XiB - (l5.2.1()ll gives the lemma in the case f^j ^ f^j . 
The proof in the case fgj > Jqj follows by swapping + and — in the definition 
ofTand.^'*. □ 

5.3 Hausdorff stability of random-cluster boundaries 

Consider the context of Proposition [F. 1.11 E^(6i) > £^(62)- Under /^^'^ 
the plus phase is dominant so the minus boundary does not affect the bulk of 
the domain. In this section, we will show that, in a random-cluster sense, the 
d~ h. boundary-cluster is small. 

Proposition 5.3.1. Let Ehs > 0. Suppose that E^(6i) > £^(62(1 - ^hs))- There 
is a constant Chs = Chs(£hs) > 0, independent of B^^^, such that with uniformly 
high Qo -probability 

^U+,~),h^Q-^ o We(6i, 62(1 - £hs))) < exp(-chs62//i). 

The proof develops the technique of truncation used in |5]. The differences in 
geometry, and the presence of a magnetic field, pose extra challenges. 

Proof of Proposition \5. 3.11 Let w denote the minimum Loo-distance between 
9yVe(l) and the origin, 

w — inf ||x||oo- (5.3.2) 

x6aWe(l) 

Let R = [ehsb2Nw/{8K)\ and let S = [(&2 - bi)Nw/{2K)\ . Define mesoscopic 
layers 

«--»('^-^-'-^). - 

The layers divide the annulus W6)(foi, 62) into S layers. Let Mj^k — ^f=j^i- The 
annulus We(62(l ~ £hs),^2) corresponds to the first 4i? layers, Mi^4ii, with Hi 
the outermost layer. The layers are essentially d — 1 dimensional, resembling 
scalar multiples of dWe; Figure [5] shows the case d — 2. The layers have been 
constructed so that: 

(i) The mesoscopic boxes in Mi form a surface separating E[;_i from HIj+i. 
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(ii) Each H/ is between 2 and 2d mesoscopic boxes thick. 

Define mesoscopic phase labels ^E* = {'^{i) : Bi<-(j) C Hi. 5) according to (15.2.11) . 
We will show that with high probability, a surface of +1 boxes separates Hi 
from H4jf. 

Given a phase label 'ip : {i : Bif(z) C Hi ^} — > {±1,0}, define the profile / 
of "0 as follows. For s G {—1, 0, 1}, let count the number of s-boxes in H;, 

/f - #{z : MkH) C H; and 4,(1) = s}. 



Let / = ifiYi=i ^'V- ^^^^ write J-{f) to denote the set of phase labels 
ip compatible with /. The number of configurations of the phase labels ip 
compatible with / is limited by the definition of the coarse graining. Surfaces 
of 0-boxes must separate the plus-boxes from the minus-boxes. The surfaces 
of 0-boxes cannot separate H; into more than /° + 1 connected components. 
Therefore the number of ways of assigning the labels in layer I is bounded by 

(^'^'1^ exp([/f -I- 1] log 2) = exp(/°0(logiV)). (5.3.3) 

We will say that a profile / is spanning if /^^^ -f /° > for Z — 1, . . . , 4i?. Recall 
from (pXTj) that K « Ari/(2d)^ ^j^g number of spanning profiles is less than 

0(N/K) 



1=1 ^ ^ 



which grows more slowly than exp(cA'^) for every positive constant c. It is 
therefore sufficient to find a constant Chs — Chs(£hs) > and a Qe-event J such 
that 

^ C {j : V/ spanning, </.;('(+'-^'''(* G F{f)) < exp(-2chs62Af)} (5.3.4) 

and J has uniformly high Qg-probability. The event J is defined as follows. 
Let J = for h ^ ho (for some ho > 0) so we can assume that h is arbitrarily 
small. We will appeal below to Proposition 15.1.11 and Lemmas 12.4.41 [5".2.2I and 
15.2.31 For h < ho, let denote the intersection of the associated Qg-events. 
Let a, d > 0. Consider three constraints on the label profile: 

s 

Y^f^ ^N^-'/Vk, (5.3.5) 



1=1 



max f-^ <^a(—] . (5.3.7) 



For J £ J we will check the inequality in (|5.3.4p in four parts. We will show: 

(I) For any Chs > 0, if (|5.3.5p fails then the inequality in (I5.3.4p holds. 

(II) For any a > 0, if (j5.3.5p holds but (15.3.61) fails, then the inequality in 
(|5.3.4p holds if Chs is sufficiently small. 
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(III) For any a > 0, if and ([O:^ holds but fails, then the 
inequality in (|5.3.4p holds provided that a and Chs are sufficiently small. 

(IV) If (|5.3.5p - (|5.3.7l) hold, then the inequality in (|5.3.4p holds if a, a and Chs 
arc sufficiently small. 

For part (I), choose / such that (15.3.51) fails. Inequality (|5.3.5p is an upper bound 
on the volume of bad boxes. We will apply Lemma 12.4.41 for ip e J-{f) with 
Ecg = 1 and n — f? > N"^^^ j \fK. The positive terms in the exponential 

in Lemma 12.4.41 are 

/3|a±A| +/3/i|A| = 0{hiN'^-^). 

The absolute value of the negative term is greater than CcgK{N'^^^ j \fK) which 
is a higher order of TV p.l.ip . For h sufficiently small, 

= V-) exp(-CegnK/2), 

and so by (|5.3.3p . 

e ^(/)) ^ exp(nO(logA^) - c.^nKj^) . 

Whatever the value of the constant Chs > 0, for h sufficiently small the right- 
hand side above is less than exp(— 2chs&2^)- 

Now to part (II). Inequality (j5.3.6p is an upper bound on the volume of 
minus phase. We can expect the majority of boxes to have label +1 because of 
Proposition 15.1.11 

Recall the symbol Ecg used in the definition of the coarse graining. Let V' 
denote a label configuration. Let n+ count the number of phase labels = +1 
such that IB /<■(«) is also £cg-good. Similarly for n~ . Let count the number of 
£cg-bad boxes in Hi^g. Note that 



1=1 

Let M(estb) refer to the ^ ''-event in Proposition 15 . 1 . Il Under M(estb), 

J2 >hi-hi- 2e,nM + 0{K-'). (5.3.8) 

For each box Bi<-(z) counted by n~ there are at least (1 — £cg)K'^m* vertices 
in B^(z) n Bi<-(i). For each box B/^ (z) counted by n+ there are at most (1 + 
ecg)K'^m* vertices in nB/f(i). There are at most nPK'^ vertices in Scg- 

bad boxes. The remaining vertices lie in clusters with diameter less than K/2. 
Small clusters only interact weakly with the magnetic field. If an open cluster 
has volume of less than {K/2)''' then the odds of it taking plus spin are at 
most ex.p{l3h{K/2y) to 1 (|2.2.ip . Conditional on * = V, with probability 
1 - exp{-b2N), 

^— ^ <j{x) ^n+{l + £,g) + ^ - n- (1 - seg) + (^) ' 0{K-'). (5.3.9) 
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Figure 2: Part (III): The shaded region U corresponds to the blocks with phase 
label —1. The arrows indicate paths from —1 phase labels to +1 phase labels. 
The black lines marking the intersection of the shaded region with M„i and H/ 
correspond to Ai and A2 in Lemma [5.3.12[ respectively. The intersection of dU 
with Hi^rn corresponds to the set S. 



By the argument from part (I) we can assume that 71'' ^ N"^ ^ jyfK. The 
number of £cg-good boxes in A is therefore 



By (io:a - (io:TU)) . 



n++n- = (^) [6^-6f + 0(X-i)]. (5.3.10) 



n- ^ (f) [£stb6^ + ecg(&^-fef)/2 + 0(i^-i)]. 

Taking e^^h = a/2 and Ecg sufhciently small, we see that (|5.3.6p holds with high 
probability. Taking 2chs < inin{l, Cstb(£stb)} we have completed part (II) of the 
proof of (|5Xil) . 

Now for part (III). Choose / such that (|5.3.5p - (j5.3.6p are satisfied but (|5.3.7p 
is not. Let \ ^k<l<m<n^'^R and let 1/' G ^{f)- See Figure [2J 

We can apply Lemma 15.2.21 with F equal to the set of mesoscopic boxes in 
H/c+i^„_i in the neighborhood of a or a —1 box, 

F-IJBk(i), (5.3.11) 

/ = {^ : 3j, BK(^),BA'(i) C Hfc+i,„_i, ^(j) < 1 and \\i - ^ !}• 

Note that F is not necessarily connected. 

Let V' G •^(/)- Recall that the quantities /af, /a/, /g}, //^ in the statement 
of Lemma [5.2.21 count the number of +1, and —1 phase labels in 9/, and the 
maximum flow from plus to minus labels in F. To estimate the quantity we 
need a geometric lemma. Recall the definition of d p. 3. II) . 
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Lemma 5.3.12. Suppose that < ai < 02 and U C Aq. Let Ai — UC\d'We{ai), 
i = 1,2. Let S denote the portion of dU contained in the closure o/Wg (01,02). 
Then 

n''-\s) ^ d-^/^[n''-\Ai) - n''-\A2)]. 

Proof. Let P denote the linear projection x — >■ (01/02)2;. S must separate 
Ai\PA2 from {P~'^Ai)\A2. The sm-face areaof Ai\FA2 is at least n'^-\Ai)- 
'H'^^^{A2) as F is a contraction. Let x e dWe. Recall that We ^ AeH W2^; by 
the symmetry of yV27r, the angle between the vector x and the normal vector 
n^»(x) is at most cos~i(d-i/2). □ 

Lemma [5.3.121 implies that there is a constant c G (0, 1) such that 

fr>cf^'~c-\f,''+f°). 

Given d, if a is sufficiently small then we can choose k < I < m < n such that 

fk+.fk' ^ {b2N/Kf-^ X dCcutc/(6Ceut), 

fi + fr' ^ {h2N/Kf-' X dcVe, 

f-^ ^ {b2N/Kf-' X d, (5.3.13) 

/n +/n ' < ib2N/Kf-^ X dCcutc/(6Ceut). 

Lemma 15.2.21 gives 

ip'/+-^^\lnt^, n Ext^) ^ exp (Ceut(A" + f^' + fn + .fn')K'''' - CcutfrK'^') 
^ exp(-dccutc(fe2A^)''"V2). (5.3.14) 

We complete part (III) by checking that the right hand side of (|5.3.14p . when 
multiplied by the size of the set J-{f) [cf. (|5.3.3p and (|5.3.5I) ] is less than 

exp(-2chs&2A^)- 

Before we start part (IV) of the proof of (j5.3.4p , we will give an isoperimetric 
inequality for dWe ■ The surface dWe is d — 1 dimensional, so subsets of dWe 
have d — 2 dimensional boundaries. 

Lemma 5.3.15. There is a positive constant u — u{6) such that for any A C 
dWe, 

C^-^idWe] ^2=^^ [dA]^u{C [A]) 

Proof. First consider the case 9 ~ 27r. Let w denote the minimum Loo-distance 
between dWg and the origin (j5.3.2p . By convexity, dWe lies inside the L2- 
annulus with inner radius w and outer radius wd. 

Consider the projection P of dWg onto the unit sphere S"^"^. Associated 
with P are two Radon-Nikodym derivatives, one for the d — 1 dimensional 
Lebesgue measures on the domain and codomain, and one for the d — 2 di- 
mensional Hausdorff measures on the domain and codomain. By symmetry 
and convexity (see [TH Theorem 2.2.4]) both Radon-Nikodym derivatives are 
bounded away from and 00. The result follows from Levy's isoperimetric 
inequality for the unit sphere. 

A similar argument works when < 9 < n. Consider a projection from We 
to the {d — l)-dimensional unit ball. □ 



31 



Now for part (IV) of the proof of (|5.3.4p . Let / denote a spanning profile and 
let ip S J^if)- Choose k and n such that R ^ k ^ 2R ^ n ^ 3R. We will apply 
Lemma [5 . 2 . 31 with F defined according to (|5.3.1ip . 

If (j5.3.7p holds with a sufficiently small then Lemma [5.3.151 implies that for 
some constant c > 0, 

n-l 
l=k+l 

Let 

Lemma 15.2.31 gives 

n Ext^) (5.3.16) 
< exp ( Ccut(A° + /,7^ + f^ + fn')K'-' + f3h\r\4'' - d] • 

V i=k+i ) 

The term corresponding to the magnetic field in (|5.3.16p is bounded by (|5.3.7p . 
If a is sufficiently small, 

;=fc+i ;=fc+i 

With reference to (|5.3.3I) and the assumption that h is small, the number of 
ways of choosing (■)/'(«) : IB ;<-(«) C Mk,n) is at most 



X^exp(/°0(log A^)) < exp ( (/° + /°)0(logiV) + J ^ 
Thus if h is sufficiently small 



n-l N 

m 



4 

l=k \ l=k+l 



+ e Hf)) ^ exp ( 2Ceut(A" + fk' + fn + fn')K''-' - ^ ^ ^ ■ 

\ l=k+l / 

(5.3.17) 

In our notation, O (4.40)] states that if / satisfies (|5.3.5p - (|5.3.6p with a suffi- 
ciently small then there is a positive constant c such that 



mm 

R<k<2R 



mm 

2i?,< 



\2C,,,t{fk+fk')K''-'-\ 9i\ ^ -cb2N, and (5.3.18) 
I i=k+i J 

n<3R. Z 

y i=2R+i J 



Part (IV) of the proof of inequality (15.3.41) follows from (|5.3.17p and (|5.3.18p . 
This completes the proof of Proposition 15 .3 . Il □ 
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We will give three analogous results below. They follow, mutatis mutan- 
dis, from the proof of Proposition 15 . 3. ll Consider first the context of Proposi- 
tion [5X101 

Proposition 5.3.19. Recall the event C defined by (|5.1.9p . Suppose that 6i = 
and 62(1 — £hs) > + £c- Given Sc and Ehs, with uniformly high ^0 -prohahility 

^j;^'''id-A ^ Wfl(6i, 62(1 - Ehs)) I C) < exp(-chs62//i). 

Proof. In part (II) of the proof of Proposition 15.3.11 replace Proposition 15.1.11 
with Proposition 15 . 1 . lOl □ 

In the context of Proposition 15.1.111 the plus phase is dominant and so the 
minus boundary does not affect the bulk of the domain. 

Proposition 5.3.20. There is a constant cf^^ = c{jg(ehs, S^iax) > such that 
with uniformly high Qg -probability 

^M-^+)^^9-^ ^ + Ehs, &2)) < expi-c'Uh). 

Proof. Proposition 15.3.201 differs from Proposition 15.3.11 in that it shows that 
the inner (rather than the outer) boundary condition has limited influence. Let 

, /, 2(AR-l)K , 2(4:R+l~l)K\ 

2(1 -1)K 21K\ 
+ + ^-4i.+ l,...,^, 

with R = \eY,sNw/{m)\ and S = [(62 - 5i)iVu'/(2i^)J . The proof of Propo- 
sition 15.3.11 mutandis mutandis, shows that a surface of +1 boxes separates 
Hi from MiR. In ([53:^ - ((5X7)) . replace h2N/K with N/K. Proposition |5JJJJ 
replaces Proposition 15 . 1 . ll in part (II) of the proof. 

Notice that the proof of part (III) has become slightly more flexible; the 
additional flexibility will be important below in the proof of Proposition 15. 3. 2T] 
The quantity cj^^ is allowed to depend on iJ^ax- This means that Lemma [5.2.31 
can be used in place of Lemma 15.2.21 the extra term due to the magnetic field 
can be controlled by taking a, and therefore |r|, sufficiently small. 

□ 

Consider the context of Proposition 15.1.81 The minus phase is dominant so 
the plus boundary does not affect the bulk of the domain. 

Proposition 5.3.21. Suppose that E^(6i -I- ehs) < £^(62)- There is a constant 
cfj'g = c[jg(ehs) > such that with uniformly high -probability 

^M+,-)^(^9+^ o W9(6i +£hs,&2)) ^ exp{-cUh). 

Proof. The proof can be obtained from the proof of Proposition 15 .3 . 2U1 by swap- 
ping the roles of plus and minus. The sign of the magnetic field has to stay the 
same, but, for example, in (15.3.611 replace fj^^ with /;"^^, etc. Proposition 15 . 1 .81 
replaces Proposition 15 . 1 . iT] in part (II) of the proof. □ 
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5.4 HausdorfF stability implies spatial mixing 

In Proposition l5.3.2T] we showed that under Qg[fi'^^^' with high probabihty, 
a surface of —1 mesoscopic blocks separates the region Wg(&i + ehs,&2) from 
the plus boundary d'^A. In two dimensions, by planar duality, this implies 
that there are no Ising spin-clusters connecting (9+A to Wg(6i + ehs,&2)- By 
the Ising model's domain Markov property, and monotonicity, we can compare 

Q^Ma^'"'"''"] to Qe[//~ in W^(6i + eh„ 62). 

In contrast in higher dimensions, especially when close to the critical tem- 
perature, the Ising spin-cluster associated with d~A under /i^'^^ may be 
much larger than the cluster associated with d~A under the random-cluster 
representation ' cannot make such a comparison. Instead we will 

appeal to a spatial mixing property, stated below as Proposition 15.4.11 

Much is known about the spatial mixing properties of the Ising model in 
the absence of a magnetic field. The difficulty here is that the magnetic field is 
acting to weaken the dominant phase. 

We conjecture that Proposition l5 . 4.11 and the coarse graining property, holds 
for all (3 > Pc- If that is the case then Theorem 11.2.31 holds up to the critical 
point. For simplicity we will reuse the constants Chs, c(jg and c((g, adjusting their 
values if necessary. 

Proposition 5.4.1. There is a finite /3q such that if (3 > /3q then for Ehs > 0, 
with uniformly high Qg -probability: 

(i) ForxeWe{bi,b2{l-2ehs)), 

^■U+,-),h = i\d-A^ We(&i, 62(1 - £hs))) 
^ >f'/+'+^'\a{x) = 1) - exp(-chs&2//i). 

(ii) ForxeW'g{bi+2ehs,b2), 

^^^^^ = i\d-A^ Wg{bi + eh., 62)) 

(ill) If E«(6i + 2eh3) < £%b2), then for x e Wg{bi + 2£hs, &2), 

^JX+-),h ^^(^^ = i\d+A^ Weibi + £hs, 62)) 
^^'/-^-^^'\aix)^l)+eM~c'Uh). 

(iv) If 61 > Bl for X e W;,(fei + ehs&2, 62), 
^'^^-^^\a{x) = 1) 



M-i^AC" ^ {BlY +exp(-chs62//^). 

We (61) / 



The statement of Proposition 15.4.11 is fine tuned to suit our needs — we have 
only considered spatial mixing in Wulff-shaped regions. Also, the restriction in 
part (iii) is stricter than necessary. 
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We will prove Proposition 15.4.11 after first showing how it can be used with 
the results of Section 15.31 By part (i) , in the context of Proposition 15.3.11 with 
uniformly high Qe-probability, for x S Wg{bi,b2{l — 2ehs)), 

\^i■'/+-'>'\aix) = 1) - ^ij:+''\a{x) - 1)1 s; 2cM-Chsb2/h). (5.4.2) 

By part (ii), in the context of Proposition 15.3.201 with uniformly high Qg- 
probability, for x £ Wg(6i + 2£hs, ^2), 

\//-'+^^\aix) = 1) - ^ii+^'\aix) = 1)1 2cM-c{Jh). (5.4.3) 

By part (iii), in the context of Proposition 15.3.^ with uniformly high Qg- 
probability, for x £ Wg(6i + 2£hs, ^2), 

l/Z-^'-^-^aix) = 1) - = 1)1 ^ 2cM-cVh). (5.4.4) 

Suppose that 62 > 61 > -Bfoot ^^'^ consider x £ Wg{bi + £11362, &2)- By part 
(iv), and by Proposition 15.1.11 applied to 'Wg{bi) with £stb = B^/^i, with 
uniformly high Qg-probability, 

^ exp(-Cstb6i//i'^"^) +cxp(-chs62//i)- 

Proposition 15.4. II is also relevant to the conditioned measure defined in (I5.1.9|) . 
With bi — 0, the total variation distance between fi/^ ' and fJ-/^~' ( ■ | d~ A 
We (62(1 ~ ^hs))) is bounded by Proposition 15.1.1^ monotonicity, and Proposi- 
tion [53?lH Thus by part (i) of Proposition l5.4.11 for some constant c > 0, with 
uniformly high Qg-probability for x £ Wg(62(l ^ 2£hs)), 

\(ii-^\a(x) - 1) - MA+''(a(a;) - 1)| exp(-c//i). (5.4.6) 

If b2 > bi + 2ec then part (iv) can be used to compare /i^''^' and 

'"w~(0b2)" ^ ^ Wg(6i +£hs&2,62), with uniformly high Qe-probability for 
some positive constant c > 0, 

l/xf +-)'^a(x) = 1) - ft^$^^^_-^{<T{x) = 1)1 ^ exp(-c//i). (5.4.7) 

Proof of Proposition \5.4-l\ We will prove part (i); the other parts are similar. 
By monotonicity it is sufficient to show that for x £ Wg{bi, 62(1 — 2£hs)), 

^'WelblM{l-e^.)) (^(^) = 1) ^ ^'^itub2{l-e^.)) ('^(^) = 1) ~ exp(-Chs&2//l). 

Let A denote the event that the inner- and outer-boundaries of We(62(l ^ 
2£hs), ^2(1 — £hs)) are separated by a set of plus spins blocking all paths between 
the two. Of course, the set of plus spins only needs to block paths composed 
entirely of edges with J(e) = 1. By monotonicity 

/^S[!il-e„.))(^(^) - 1 I ^) 55 MwllbitMl-e.,))^^^^) = 1) 

so we need to show that 

A^wlti2(i-e,,.))(^) > l-exp(-chs62//i). 
We will do this using a stronger coarse-graining property. 
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Definition 5.4.8. Consider a boxMKii) C A. // 



(i) Bx(j) is Ecg-good and 

(ii) the a-spin clusters composed of vertices with spin — cr(Bj^(i)) intersecting 
Bi<-(i) have diameter at most K/2 

then say that Bx(i) is £cg-Ising-good. 

With p G (Pc,l) fixed, as /3 — >■ cx) the annealed random-cluster measure 
Q[(/)'^'''] converges weakly to product measure with density p; the density of 
edges with J(e) 1 but ^(e) = goes to zero. Taking Kq large, and then 
taking /? large, we can make the Q[(p^''*''"]-probability that MKgii) C A is £cg- 
Ising good arbitrarily close to 1. By a standard renormalization argument we 
can find f3o, Ki and c > such that if /3 > /3o and K > Ki then B^f (i) C A is 
£cg-Ising-good with Q[</?a '°]-probability 1 — exp(— cif). 

The result now follows by adapting the proof of Proposition 15.3.11 Sub- 
stitute '£cg-Ising-good' for 'Scg-good' in the definition of the phase labels and, 
because of the free outer boundary conditions, use Proposition l5 . 1 . 1^ in place of 
Proposition 15 . 1 . Il If the profile / associated with the phase label configuration 
^ is not spanning then the event A holds. 

Parts (ii) and (iii) of Proposition 15.4.11 follow from the proofs of Proposi- 
tions 15.3.201 and 15.3.211 respectively, by substituting 'Scg-Ising-good' for 'Scg- 
good'. Part (iv) follows from the proof of Proposition 15 .3 . l"9l with high proba- 
bility there is a surface of plus spins separating the inner- and outer-boundaries 
of WoibiM + under ■ \ J^^^^^^ M], dC^ ^ {B^). □ 

5.5 Spectral gap of the dynamics 

The Glauber dynamics for fJ.'j^'''^ can be studied by introducing a block dynam- 
ics. With ebiock > 0, let n = [(62 — fei)/ebiockJ — 1- Consider a sequence of 
overlapping annuli that cover A, 

Aj = Wg{bi + (j - l)£biock, bi + {j + l)ebiock), j = 1, 2, . . . , n - 1, 
A„=W;,(6i + (n-l)ebiock,62). 

Consider a block dynamics for /i^'"'''* with blocks Ai, . . . , A„; update each block 
Aj at rate 1, resampling the block conditional on the configuration restricted 
to A\Aj. 

Lemma 5.5.1. Fore > 0, i/ebiock andho are sufficiently small andO < h < Kq, 

gap(A, C, h) ^ exp(-e//i'^-i)gap(A, {Ai, . . . , A„}, (, h). 

Proof. We is a subset of yV27r, so we can assume — 211 without loss of generality. 
Let yl, . . . , j'l'^ I denote an ordering of the vertices in Aj such that the angle 

between yf and ei is increasing with i. For each vertex yj in block Aj consider 
the edge-boundary between {yJ, . . . ,yf} and {yf^i, ■ ■ ■ , J/j'^ |}! ^ denote the 
maximum (over i — 1, . . . , | Aj | and j = 1, . . . ,n) size of the boundary. Given 
B^^^, Lh'^-'^ = O(ebiock) as Sbiock 0. 
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As noted in [TB^, the proof of [HI Theorem 2.1] imphes that for some C, c > 0, 
gap(A, C, h) ^ ^^^^^l^^g&p{A, {Ai, . . . , A„}, C, h). 

Choose Ebiock so that CL < e/h'^~^. □ 

We are now in a position to extend [H] Propositions 3.5.1-3.5.3] from the 
Ising model on to the dilute Ising model on Z'' with d ^ 2. 

Proposition 5.5.2. Let bi ~ and e > 0. With uniformly high Qg -probability, 

gap(A,+,/i) > exp{-e/h'^^^). 

Proof. Let {(7t)t^o denote a copy of the block dynamics Markov chain. The 
graphical construction can be extended to the block dynamics by coupling from 
the past: if block Aj is to be updated at time t, use a copy of the regular graph- 
ical construction in Aj over the time interval (— oo, 0] with boundary conditions 
(Tt_ to produce the new configuration at- By monotonicity, at is an increasing 
function of the initial configuration cto. 

When t is sufficiently large, at is independent of ctq; at then corresponds to a 
sample from the equilibrium distribution Let {at^"^)t ^ o denote a copy 

of the block dynamics Markov chain started in equilibrium. 

We will show that with uniformly high Qg-probability, the probability that 
ai = crj''™ is bounded away from zero. This implies that the spectral gap of 
the block dynamics is bounded away from zero and so the result follows by 
Lemma 15.5.11 

Say that an update of block Aj at time t is good if the update maps all 
configurations that agree with ctj^'" on A\ U^^^^A^ to configurations that agree 
with CTj'^™ on the larger set A\ U^~-J A^ = Wg{bi + jEbiock, &2)- By monotonicity, 
if at- agrees with cr"^™ ~ Ma'^' on A \ U^^;^ A^ then 

^w;,(f)i+0-i)ebiook,&2) ^''t ^* '^w;,(6i+(i-i)ebiock,f'2)- ya.0.0) 

Inequality (j5.4.3p used with the sandwich (j5.5.3p gives a lower bound on the 
probability that the update is good; at and a't'^™ agree on Wg(bi + j£biock, ^2) 
with probability at least 1 — 2|A| exp{—c'^^/h). 

With probability exp(— n)/n! there is an uninterrupted sequence of updates 
on An, A„_i, . . . , Ai in the time interval [0, 1]. If all the updates are good, which 
occurs with probability at least 1 — 2n|A| cxp(— cj^^A^), then ai = a'l'^™. Note 
that n ^ _B^-jjjx/£biock so the probability of seeing such a sequence of updates is 
bounded away from zero uniformly over 62 € [0, S^iaxl- ^ 

Proposition 5.5.4. Let hi ~ and e > 0. One can choose 62 slightly larger 
than such that with high Q^g -probability, 

gap(A,-,/i)^ exp(-£//i''"i). 

Proof. Take ebiock according to Lemma 15.5.11 Taking n = [i^c/^biockj — 1, 
choose 62 G {BI, {n + 2)ebiock) such that E^((n — l)ebiock) < £^(62)- We can 
then follow the proof of Proposition 15.5.21 The boundary conditions in (j5.5.3p 
should be changed to (— , — ) on the left-side and (+, — ) on the right-side. The 
probability of a block update being good is then bounded below using inequality 
(ICTl) in place of (ICT^ □ 
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Proposition 5.5.5. Consider the case bi > B^. Let e > 0. With uniformly 
high -probability, 



Proof. Say that an update of block Aj at time t is good if the update maps ah 
configurations that agree with CTj^™ on We(6i + (j — l)£biock) = A \ U"^^ to 
configurations that agree with ctj'^™ on A\ U"^j_|_]^Ai. The boundary conditions 
in (|5.5.3p should be changed to (+, — ) on the left-side and (+, +) on the right- 
side. Inequality (15.4.21) shows that updates are good with high probability. If 
there is an uninterrupted sequence of good updates in the order Ai, . . . , A„ in 
the time interval [0, 1] then cri = a'^^™ . □ 

6 Space-time cones and rescaling 

In this section we will turn the heuristic description of plus-cluster nucleation 
from Section 11.31 into a proof of Theorem 11.2.31 We will apply the results in 
Section [S] with two values of 9. 

We will take 6* e (0, tt) to denote the argument of in the statement of the 
theorem. We will consider regions with the shape Wg{b) for b e [B^^^, B^^^^.]. 
The lower bound Bf-^^^ will be chosen to maximize the rate of nucleation of plus 
clusters. We will take -B^ax to be the minimum value such that a translation of 
W2^(1.01S2'^) fits inside We(S^ax); see parts 1 and 2 of FigureH By Proposi- 
tion HXl W27r(5c'') lias diameter of order 1 as /3 -> oo and 6* 0; yVeiBf^^^) 
must have volume of order d~^. The Wulff shape yVe(fo) has volume 6*^ so iJ^^^x 
must be of order Q^'^l'^ . By (14. 3. 21) the probability of the event conditioned on 
in the definition of Qe is 



This gives the density of nucleation sites in . We will show that at these 
nucleation sites, droplets of plus phase form at the rate exp(— E^//i''^^). 

We must then show that the clusters of plus-phase can spread out from the 
sheltered nucleation sites. We do this by considering the full Wulff shape Y^2t!- 
In areas of typical dilution, sufficiently large Wulff-shaped droplets of plus phase 
expand with high probability. With reference to (j5.4.2p we will take B^^^^ large 
so that fJ'^~'^g2,T ) provides a good approximation to the equilibrium measure 
/i'^''* in a neighborhood of the origin. 

6.1 The graphical construction in space-time regions 

Before we give the proof of Theorem ll.2.3[ we need to extend the Ising dynamics 
to allow the size of the graph to change with time. With Fq, Fi, . . . , r„ C Z''' 
and to < ti < ■ ■ ■ < tfi+i, consider the space-time region 



gap(A,(+,-),/i) exp(-e//^'^-i). 




(6.0.6) 



n 




(6.1.1) 



i=0 



The graphical construction for the Ising model /^^' 
can be extended to F. 



described in Section [^31 
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Figure 3: Illustration of a Q-catalyst. The nucleation event occurs at rate 
exp(— E^//i''~^) and consists of the following steps. (1) A droplet of plus phase 
with the shape W0{B^^^) forms in a region of high dilution that resembles 
Wg (S^ax) [Proposition l6.2.2] . (2) The droplet expands in the sheltered region to 
cover a copy of W27r(l-01-Bc'^) [Proposition l6.3.1j . (3) The droplet of plus phase 
spreads to the right [Proposition 16.4.2] and (4) expands to cover W27r(^in'in) 
[Proposition 16.3.11 with 9 taken to be 27r] . There is now a droplet of plus phase 
at the center of a Q-conductive site of the rescaled lattice. The droplet ex- 
pands (5) to cover W27r(i3^^x) [Proposition 16.5.1] which contains neighboring 
Q-conductive sites. 
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(i) Let s denote the start time. 

(ii) Let ^ denote an initial configuration compatible with boundary conditions 
C at time s, i.e. if s G [ti,ti+i) then ^ G Ep . 

(iii) Let crp|_^.^ = ^. 

(iv) If a vertex x is added to the dynamics at time ti (i.e. a; G \ r^-i) then 
the spin crp| ^.j. (x) is taken to be £^{x) to match the boundary conditions. 
The spin at x may then change with each arrival of the corresponding 
Poisson process. 

(v) If X is removed from the dynamics at time ti (i.e. x G Ti^i \ Ti) then 
the spin at x is immediately switched to ^(a;) to conform to the boundary 
conditions. 

The graphical construction of Pj allows us to link together the Ising dynamics 
run in overlapping space-time regions. This can be used to chain together the 
different steps involved in the growth of a region of plus-phase. 

Remark 6.1.2. Consider two space-time regions such that the top layer of the 
first region covers the start of the second region: 

F = ST(Fo, Fi, . . . , F„i; tg < <i < • • • < t,n+i), 
A = ST(Ao, Ai, . . . , A„; mq < ui < • • • < m„+i), 
Tm = Ao and uq = < t,n+i ^ "i- 

"ruA,-,/l;M„ + i ~ "A,-,h;n„ + r 

6.2 Droplet creation in a Summertop cone 

Let e e (0, tt) and (5 > 0. By Proposition lSXil we can choose B^^^^ e {B^,B^^^^) 
such that with high Qe-probability, 

gap(We«i„), h) ^ exp(-<5/(2/i^-i)). (6.2.1) 

Let A = Wfl(i?J^ji,j). Heuristically, we expect critical droplets to form in A at 
rate exp(— E^/ft,''^^). Let jl'j^ '^ = Ma ''*( ' I ^) denote the conditional measure 
defined by (15X91) with ec = {B^,,;^ - B^^)/3. 

Proposition 6.2.2. With high Qg -probability we can construct a random vari- 
able (7 ^ fi'j^ such that the event {a = '''a - /i-exp((5/h''-i) J' probability 
exp(— E^/ft,''^^) and is independent of the value of a. 

Proof. Taking a = in the last inequality in the proof of Proposition 15.1.101 
we can assume that //^'^''*(C) ^ 2exp(— E^/Zi''"^). By (j6.2.ip and a Markov 
chain mixing inequality (i.e. |15l (59)]) the total variation distance between 

'^A^-MoMS/h"-^) ^^^^ ^A~''' ^^'^^ ^^^^ exp(-E«//i''-i). □ 
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6.3 Growing in a Summertop cone 

In this section we will use the "inverted space-time pyramids" of |16) to show 
that under Qg, droplets of plus phase tends to expand from WglBf-^^^^) to 
We(Bf„^^) with high probability. 

With (5 > 0, and with reference to (|4.3.3I) and (|6.1.ip . consider the space-time 
region 

V = V(i3,^^in, i?niax) ^7 ^) := ST(Ao, . . . , A„; to < ■ ■■ < tn+l), 

A, a;+I^^(^-)I, U ■.^^eMS/h'-'), n := |W,(i?Ljl " |W«(i?Lj|. 

For each i let ^a'^''' = A'aT''*! ' I with ec = (-B^^ - Bl)/i. 
For rj G E^^^^e ^ consider the event 

^ /^O.'' — ^*"'+ \ 

For r] ^C, Gn describes the plus phase spreading from We(i?j^j;jJ to ^^eiB^^y) 
in time The event Gj, depends only on the elements of the graphical con- 
struction contained in V. Here is an extension of (THl Proposition 3.2.2] to the 
dilute Ising model. 

Proposition 6.3.1. There are positive constants (5o,C, Gi such that if < 
S ^ Sq then with high Qg -probability 

I Fj{G,) dfii:-\Tj) ^ 1 - CcM-Ci/h). 

Proof of Proposition \6.3.1[ We will assume that J belongs to a certain event 
with high Qg-probability; the set is defined implicitly by our use of results from 
Section [5l 

For ^ = 1, . . . , n and C e S^^^ let G^ = {^^vT-Vm^+i = 4V-,?>;t,+i }: 

n (^p GV^ =^ G,. 
By monotonicity G^ C G\. It is sufficient to show that for each i, 

J Pj{Gl) dAt:f (C) > 1 - Gcxp{-Ci/h). (6.3.2) 

With reference to Proposition 15 . 1 . lUl if we start the dynamics with initial dis- 
tribution /tA."'/' we expect to stay inside C for a long time. Let (o'^'^_ ^.Jf ^ ^ 
denote the Markov chain obtained from the graphical construction by suppress- 
ing any jumps from C to C^. By introducing a stopping time 

= inf{i ^ ti : crvV-Ji;* ^ '.- .h-A ^ 

we will see that the modified dynamics are likely to agree with the regular 
dynamics over the interval [ti,ti^i]. 
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Let denote the configuration obtained from a by fiipping tlie spin at x, 
and let 

9C = {(TeC:3x,CT'^eC'}. 

Proposition IS . 1 . lUl gives an upper bound on /i,^'. '''(9C). Given iJ^j^ we can find 
(5o > such that 

fi-j^-^'^idC) s; exp{-3So/h'^''^). (6.3.3) 

If the starting state C is sampled from Aa ~''' then the process ((T^ 'i ^.t)te[ti,ti+i] 
is stationary. By (|6.3.3p (cf. [16, (2.12)]), if 5 ^ (5o and h is sufficiently small, 

J Pj(Tf ^ i,+i) dAA^^'HC) < exp(-<5o//i'*-i). (6.3.4) 

For some x, = Ai_i U {a;}. We will need a bound on the effect of adding 
this extra vertex has on the conditional Ising measures (/i ]('"''').;. By the Ising 
model's finite-energy property, for any ho > 0, 

a:= inf inf inf inf //•\''*(ct(0) = s) > 0. (6.3.5) 
By the Ising model's Markov property, for ^ e C n , 

Air''(o/Atf (c) = MAr''(^(^) = -1 1 c) ^ a. 

Therefore (cf. P; (3.28)]), 

y'pj((G^)^)dAt:f(C) (6.3.6) 

=j^.ji^:-iu,. + 4'-,.;t.,,)dAtf (c) 

^ /ip^(<'i,M,.. ^<-,/.*,+JdAt:f(c) + /p,7(r;Li «;i.)dAi:f(c) 

a-iy P,7(a;v5,„,,,^, ^ 4'+,^,^^JdA^r'"(C) +exp(-5o//i''-'). 

Set 

^ = (l-7)Sc'+7Smin, 7e[0,l]. (6.3.7) 

Thus = bo < bi/3 < &2/3 < bi = Choose 6 such that W'g{b) = A^. 

By monotonicity and the invariance of the modified dynamics with respect to 
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^W^(b)' 



!yew;,(6) 



< E {^•^(^W^"(6), + >;cxp(5//i<'-i)(y) = 1) - Aw^(i(f^(y) = 1)} 
y6We(62/3) 

+ E {lP^('^wj6i/3,6),(+,-),'»;oxp(5A''-i)(2^) = 1) - Aw^(6)(^(y) = 1)} 
y6W^(62/3.f>) 

+ exp(-<5o//i''"'). 

For y G We(52/3), by Proposition 15 .5 . 2l and Markov chain mixing [TSl (59)], 

\^ji^wlib), + M:oMS/h^-')iy) = 1) - ^ 

sc: exp [-exp(V/i''-')gap(W^(5),+,/i)] /^^^(^'^('t - +) 
exp [-exp((5/(2/i'^-i))] . 
Similarly for y G Wg(62/3,fe), by Proposition 15.5.51 

l^"'('''w^(bi/3,b),( + ,-),/i;Gxp(<5//irf-i)(y) = -*-) ^ A^W^, (hi/Lb) "'"•'I 

^ exp[-exp(<5//i''-i)gap(W;,(5i/3,6),(+,-),/i)] 741^3,6) ('^ = +) 
^ exp [-exp((5/(2/i''-i))] . 
By the above 



^ E 4^5(^(2^) = +l)-'Xw(^(^) = +l) (6.3.8) 

yGWe(62/3) 

+ E /^wllbTii (^(y) = +1) - Aw^(6) (^(y) = +1) 

yew;, (62/3,6) 

+ eM-So/h''-^) + |W'e(&)| exp [- exp(J/(2/i''-i))] . 
By (jCT^ and ([5X71) . 

yeWe(62/3) 

+ E ^wltbT/ii ('^(^^ = +1) - Aw^('6'') = +1) 

!/GW;,(b2/3,b) 

^ |W^(6)|exp(-c//.). 
Inequality (jO:^ now follows by (gSSl), and □ 
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6.4 Escaping from Summertop-cones 



In Proposition 16. 3.T] we considered space-time pyramids. Consider now "space- 
time parallelepipeds" . From now on we will write 27r in place of to make it 
clear that 9 refers to the angle of the catalyst cone. Let a denote a positive 
constant and let b = l.Oli?^'^. We can find a sequence of graphs Aq, Ai, . . . , A„ 
such that 

(i) Ao-W2.(6), 

(ii) A„ = W2^(&)+aiVei, 

(iii) Ai+i differs from A^ by adding a vertex or removing a vertex, 

(iv) for any i, for some ki e {0,aN), Ai differs from We{b) + hei by at most 
a mesoscopic layer of vertices around the boundary, and 

(v) n^O{ab'^~^/h'^), 
Let 

= 0(a, b, S) = ST(Ao, . . . , A„; to < • • • < Wi), (6.4.1) 

with ti := iex.Tp{5/h'^~'^). In Figure |31 the dotted lines indicate the area swept 
out by a space-time parallelepiped that starts inside the copy of We{B^^^). 
For rj G consider the event 

r< /^o.'; — ^*"'+ \ 

Here is an extension of Proposition 16.3.11 to 0- Let eq ~ {b — B'^'^)/3 and let 

At"'" :=^Ar''(- 1 c.) with 

[ JWe(B»-l-ec) + (fc,/Af)ei J 

Proposition 6.4.2. Let be defined according to (16.4. ID with S ^ Sq. There 
are positive constants C, Ci such that with high Q-probability 

J Pj(G,)dAt"'"(ry) ^ 1 - Cexp{-Ci/h). 

Proof. We can adapt the proof of Proposition 16.3.11 showing that (j6.3.2p holds 
when, for example. A; ~ Ai_i \ {x} for some vertex x on the boundary of 
Ai_i. Let fi.t)t > s denote the Markov chain obtained from the graphical 

construction by suppressing any jumps from Ci-i to C!{_i. The only place where 
the change is important is in inequality (j6.3.6p . Recall that the spin of vertices 
leaving are set to —1. Let fi denote the measure obtained by sampling from 

and then setting the spin at x equal to —1. Let /i' = Ma~'''(' I 
By the definition of a (|6.3.5p . 

M(C)<a-V'(C), Cesr- 
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In place of (|6.3.6p we have that 

y'p^((G^)^)dAtf(c) 

=lpA<^Viu,. ^ 4-,^.....)dAi:f (c) 

The rest of the proof follows mutatis mutandis. □ 
6.5 Growth on a rescaled lattice 

In Proposition 1533] we require Bf-^^^ > B^. If in addition Bf-^^^^ > B^^^^ then we 
get the following stronger result corresponding to 16, Proposition 3.2.1]. 

Proposition 6.5.1. Let B"^^^ > B^^^ > B^^^^ and 5 > 0. Consider V = 
^(^mini ^maxi -^tt). There are positive constants C,Ci such that with high 
Q-probability 

J Fj{Gr,)dt,-^-;;'lg^^Jrj) ^l-CeM-Ci/h). 

Moreover, Ci is a function of B^^^ and Ci — >■ oo as iJj^Jjj — > oo. 

Proof of Proposition 16.5. il Taking 9 = 27r, the proof of this proposition is very 
similar to the proof of Proposition 16.3.11 Define 6..^ = (1 — 7)5^^^^ + jB^-^ in 
place of (|6.3.7p . We can then simply replace fi^~' ' with fij( ' ■ The need for the 
modified dynamics and the stopping time has disappeared; the exp{—6o/h'^^^) 
terms can be removed from the proof. 

Let Ehs = (B^Tn - B^oot)/i»B^in)- MM, and (15X51) replace 

(|5.4.6p and (|5.4.7p . respectively. We can therefore replace the term exp{—c/h) 
with 2 exp(— Chs-Binln/^)- This yields the claim that Ci — cx) as Bf^^^ oo. □ 

Proof of Theorem ] 1. 2. 3[ Let A, Ag, Co and / refer to the corresponding quanti- 
ties in the statement of the theorem. Let i5 = (A — A2)/3 > 0. 

The idea of a droplet of plus phase growing can be formalized using Proposi- 
tion l6.5.fl With reference to Proposition l6.5.11 choose B^^^ such that Ci ^ Cq + 
1. Let / = diameter(>V27r(S^Tn + !))• With reference to ((5X^ . take B'^^^ to 
be greater than 3{B^-^ + 1) and large enough that for some constant C, 

- ^ cwfWoocM-c'o/h). 

Define a collection of overlapping translations of V = V(i3^'[jj, -B^ax' ^tt): let 

V^,i V + {INx, iT), {x, i)eIJ^ X N, 

where T denotes the time from the start of the first slice of V to the start of the 
final slice of V, 

T = m^Bl:^, Jl exp(<5//j^'-i). 
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Time- wise, the top slice of Vx,i overlaps the bottom slice of V^^.i+i- We have 
chosen , I and B"^-^^^ so that 



[W2.(S2-J + INei] n W2.(B^ln) = 0, and 
[W2.(Sf„Tn) + INe^] C W2.(i?^r,J. 

If \\x — 2/11 1 = 1, then V^^^o and V^^q do not intersect at time 0, but they then 
'invade' each other: at time T, V^.o covers Vy^i. 

Say that x e Z'' is Q- conductive if, translated by INx^ the Q-event from 
Proposition 16.5.11 holds; x is Q-conductive with high Q-probability. When h 
is small the Q-conductive vertices form a supercritical site-percolation type of 
process on Z''. Let G+^x,i denote the translation by {lNx,iT) of the Pj-event 
G+ from Proposition IG . 5 . II space-time paths of G+.^j^i-events show how clusters 
of plus phase spread out once they have formed. 

We will say a; e Z'^ is a Q-catalyst if the event conditioned on in the definition 
of Qe, translated by INx, occurs. Let D denote the density of Q-catalysts 

(EMI). 

For a Q-catalyst to be effective, the edges that do not need to be closed 
should have typical dilution. Choose k minimal such that 

We(S,Lx) n [W2.(i?,tJ + klNe,] = 0. 

With reference to Figure [U define a Pj-measurable event corresponding to the 
nucleation and escape of a plus droplet, 

NuC:r + -Ji;T+cxp(A7/i''-i) = '^vjXfcei.o) -,'i;T+cxp(d7h'i-i)} ■ 

Figure [3] illustrates how Nuc^, can be written as the concatenation of the events 
described in Propositions 16.2.211^3711 in the applications of Propositions 16. 2. 21 
iXUtake the value of d to be min{5o, (A - A^)/6}. 
We will say that a Q-catalyst x is good if 

(i) X + kei is Q-conductive, and 

(ii) Pj(Nuc,) > exp(-E«//i'^-i)/2. 

Q-catalysts are good with high Q-probability. If a; is a good Q-catalyst, let 
NuCa;^i denote Nuc^^ translated iT forward in time. 

Let ConM(a;, y) denote the Q-event that x and y are joined by a simple path 
of exactly M Q-conductive vertices, and let 

A := {x £ Z"^ : X — kei is a good Q-catalyst and ConM(a;,0)}, 

Take M maximal such that Vq.sjm finishes before time exp{X/h'^^^). By a 
Peierls argument, there is a constant c > such that with high Q-probability 
{\A\ ^ cDM'^}. Assume that \A\ ^ cDW^. 

The expected number of (x, i) £ Ax {0, 1, . . . , Af} such that Nuc^^^i occurs 

is 

exp i-Et/h"-') /2 • cDM^ ■ (M + 1) > exp (^^i^) » 1- 

We can assume that Nuc^j^i does occur for some (x, i) £ A x {0,1, . . . , M}. As 
X £ A, there is path yo,yi, . . . , yM £ Z'' of Q-conductive sites from y^) = x to 
Vm = 0. The time between and Vq.sji/ is between 2MT and 3MT. 
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Figure 4: The x marks a nucleation event NuCa^.i. The horizontal axis corre- 
sponds to a path of length M = 3 from yo = x ^ A to the origin j/m = 0. The 
vertical axis corresponds to time. The arrows indicate how the region of plus- 
phase can spread to neighboring points of the rescaled lattice [Proposition l6.5.1| . 
The black dots indicate points of the rescaled space-time lattice where G+.y.^ 
occurs. The spread of the plus phase is thus bounded below by a supercritical 
directed-percolation cluster. 



The growth of the region of plus-phase along the path uq, . . . , um corresponds 
to a directed percolation cluster on the graph {0, 1, . . . , M} x + 1, . . . , 3M}; 
see Figured] By a Peierls argument, with high probability there is a space-time 
path 

{{yj,,k):k^i,..., 3M; j, = 0, jsm = M and V/c, \jk - jk+i \ ^ 1) 
such that G+^y.^ _k occurs for /c = i 4- 1, . . . , 3M. □ 

This completes the proof of Theorem 11.2.31 part (i) . Part (ii) follows by 
Proposition l4.2.4l 
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